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WELLCOME STITUTE 
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FOREWORD 


The State Water Quality Control Board, after sponsoring in 1954 an exploratory 
investigation of all aspects of the submarine outfall disposal of domestic wastes, 
contracted for research in various specific fields where oceanographic studies were 
needed. Among these studies, a major project was an oceanographic and biological 
survey of the Southern California mainland shelf undertaken by the University of 
Southern California, Allan Hancock Foundation, during the period of 1956-61. 

Early in 1961, the Allan Hancock Foundation proposed a cooperative research 
program relative to marine disposal of wastes, going beyond the previous studies 
conducted for the State Board. The principal aim was the investigation of the fate 
of organic waste and nutrient substances in the sea and their relationship to bio- 
logical productivity. The project was undertaken during the fiscal years of 1961-62 
through 1963-64 with the largest share of financial support by Demonstration 
Grants from the U. S. Public Health Service. The State Board participated in 
support of those portions of the program which would have the most direct and 
immediate applications to water pollution control problems. Dr. Leslie A. 
Chambers, Director of the Hancock Foundation, was Project Director. The re- 
search program was conducted under the supervision of Dr. Richard B. Tibby. 

The succeeding pages present the contractor’s final project report for the State 
Board, dated September 15, 1964. Background information is given in the Intro- 
duction of the report. The final chapter, Summary and Conclusions, covers the 
four aspects of waste fields in open coastal waters which were investigated. 

Printing and distribution of the report as Publication No. 29 was authorized 
by the State Board on November 25, 1964. 

The investigations reported herein were conducted under the sponsorship of 
the State Water Quality Control Board (formerly the State Water Pollution Con- 
trol Board). The investigations and their direction were under the responsibility 
of the University of Southern California, Allan Hancock Foundation. The con- 
clusions and recommendations are, therefore, those of the research contractor, and 
do not necessarily reflect opinions or policies of the State Water Quality Control 
Board. 
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UNIVERSITY OF SOUTHERN CALIFORNIA 


ALLAN HANCOCK FOUNDATION 
UNIVERSITY PARK 
LOS ANGELES, CALIFORNIA 90007 


February 2, 1965 


Mr. Paul R. Bonderson, 

Executive Officer 

California State Water Quality Control Board 
1416 Ninth Street, Room 855 

Sacramento, California 95814 


Dear Mr. Bonderson: 

We are pleased to submit this Research Report on selected aspects of the biological, chemical, and physical 
behavior of waste fields in open coastal waters, knowledge of which can be generally applied to the solution of marine 
waste disposal problems. The investigations were conducted under Standard Agreement No. 12-17 from the Cali- 
fornia State Water Quality Control Board, and under Demonstration Grant WPD-01-63 from the United States 
Public Health Service, Division of Water Supply and Pollution Control. 

The biological portion of the program was directed mainly toward establishing the effects of waste discharge 
on primary phytoplankton productivity. One of the unique features of the biological work, and the one which makes 
the results more generally applicable, was the determination of the time-sequence rather than the spatial distribution 
of biological and associated chemical events within an open, dispersing field. The measured effects, therefore, are 
independent of local differences in current systems and other such hydrographic conditions that have, in the past, 
imposed severe limitations on the usefulness of the studies and predictions based solely on spatially static surveys 
and on empirical distance-relationships. Also, it has been shown that appropriate methods for identifying a waste 
field and for following its subsequent trajectory are not only feasible but should be more widely adopted for field 
studies on biological, chemical and other changes associated with waste discharge. 

The initial success of the dye techniques in this program soon led to a greatly expanded effort in both theoretical 
and empirical studies on the phenomenon of eddy diffusion. For the most direct purposes of the governmental agen- 
cies which have supported the program, these diffusion studies have led in some cases to confirmation in the field 
of the validity of certain theoretical equations for calculating diffusion, and in other cases to new information and 
methods, both of which will increase the confidence with which engineering predictions can be made of the dispersal 
of waste fields. 

Primary productivity and eddy diffusion also are of considerable interest in their own right as fundamental 
processes in the sea, and it is hoped that the results of this research program are of some basic, as well 
as technological, significance. 


Sincerely yours, 


ae aX al 


Richard B. Tibby, Ph.D. 
Oceanographer 
RBT:ak 
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A. DILUTION AND DISPERSION OF A WASTE FIELD IN THE SEA 


1. OBJECTIVES 


The objectives of this phase of the research activities 
for the past three years are summarized as follows: 


1. To develop equipment and devise methods to be 
applied to the quantitative study of the dilution and 
dispersion of a surface waste field in the sea. 


2. To furnish dilution data for concurrent studies on 
the fate of organic and inorganic nutrients discharged 
into the sea, and their effects on plankton productivity. 


3. To investigate some of the current statistical 
models of turbulent eddy diffusion and to substantiate 
the applicability of these models to the cases of in- 
stantaneous and continuous release of dye from a point 
source in the sea. 


4. To determine the effect of wind, water column 
stability and sea conditions on the rate of physical 
dilution of a waste and on related diffusion parameters, 
with particular attention to the role of vertical eddy 
diffusion in the overall dispersion process. 


5. To apply the results of the dye studies to the 
actual dispersion of a large scale surface waste field in 
the sea. 


2. EQUIPMENT AND METHODS 


Instrumentation 


During this investigation considerable time and ef- 
fort were expended on the improvement and modifica- 
tion of tracer detection equipment, design of an under- 
way sampling system, and the study of the physical 
and chemical characteristics of rhodamine-B dye used 
as an external tracer. The results of this phase of the 
research were reported in detail in the second Annual 
Report by the Allan Hancock Foundation (1). A brief 
review of the more important aspects of this work 
follows. 

Two G. K. Turner Model 111 fluorometers equipped 
with continuous-flow doors were used throughout this 
investigation. The fluorometers as received from the 
manufacturer were not suitable for immediate field use. 
The major problems encountered were: 


1) light leaks around the flow door and through the 
back of the instrument. 

deflection of the forward light path when the flow 
door was subjected to small normal loads. This 
deflection was due to the fact that the mounting 
bracket for the mirror assembly and cuvette 
holder was anchored to the door only on one end, 
thus normal loads on the door deflected the 
mirror assembly. 

condensation on the cuvette. On several occasions 
condensate was observed to flow down the face 
of the cuvette. The presence of condensate re- 
sulted in erroneous indicated fluorescence levels, 
especially at high instrument sensitivity settings. 


2) 


3) 


C1) 


4) the location of the sensitivity control knob inside 
the fluorometer. This arrangement is suitable for 
laboratory work with the non-flow door. However, 
when used in the field with the continuous-flow 
door arrangement, it becomes an awkward and 
time-consuming task to change sensitivity. 


Some of these problems were not serious and were 
easily solved. A light shield was built to house the 
fluorometer, thus minimizing the possibility of light 
leaks. Ventilation openings were provided for instru- 
ment temperature control. Condensation was minimized 
by forcing a continuous supply of dry air through the 
cuvette chamber by use of a small diaphragm pump 
and drying column. Deflections of the mirror assembly 
were reduced by anchoring both ends of the mounting 
bracket and securing a piece of \ inch flat steel bar 
stock diagonally across the flow door to make it more 
rigid. Both fluorometers were modified by addition of 
an external sensitivity control knob and dial located 
on the control panel above the cuvette chamber. By 
this means of external control it was possible to make 
rapid sensitivity changes without reaching inside the 
unit. This modification also extends the range of sensi- 
tivity settings to any desired value between zero and 
maximum instead of the 1, 3, 10 and 30 X values 
obtainable with the unmodified unit. Setting repeat- 
ability is provided by addition of a geared indicator 
dial attached to the control knob. 


For field use the instrument is equipped with a pipe 
inlet-outlet manifold which was designed to allow 
continuous flow through the sampling chamber or 
alternately through a by-pass line. An auxiliary inlet 
line was provided in order that standard solutions of 
known tracer concentration could be passed through 
the fluorometer while the by-pass was in operation. 
This arrangement allows spot calibration checks to be 
made at any time during the sampling process. Figure 
1 shows the modified fluorometer with light shield, air 
pump and drying column, as used in the field. 


In order to sample continuously from a moving 
reference platform, it was necessary to design and 
construct an underway sampling device. The sampler 
consisted of a 4 foot section of 2 inch galvanized pipe 
with 14 inch steel plate fins with adjustable flaps. A 
2 inch reducer attached to a ¥% inch section of pipe is 
located below the body of the sampler and the open 
end is oriented in the direction of flow. The section of 
pipe leading from the inlet passes up through the main 
body of the unit; the end of this line is beveled in order 
to fit snugly inside a % inch polyethylene sampling 
line. This design allows utilization of full impact pres- 
sure of the water medium to minimize sampling pump 
lift conditions during underway operation. 

The sampler was designed so that it could be towed 
behind a moving boat. The tail fins were adjusted to 
provide stability and keep the sampler submerged at 
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all times. A calibrated stranded wire and winch were 
used to lower or raise the sampler to the desired depth. 
A 50 foot section of polyethylene hose attached to the 
sampler was connected to a small positive displacement 
pump which in turn discharged directly into the pipe 
manifold attached to the fluorometer. Power for both 
pump and fluorometer was provided by a small 1000 
watt portable generator. 

The sampling equipment was assembled and mounted 
on an 18 ft. launch from the R.V./Velero IV, the U.S.C. 
marine research vessel. With this arrangement, it was 
possible to move personnel and equipment to any de- 
sired oceanic area. The launch was lowered for diffusion 
work, thus freeing the Velero IV for additional studies. 
The launch had the added advantage of small size and 
shallow draft, important factors when working in 
surface dye fields. 


The speed of the launch was estimated by use of a 
Walker-type speed indicator attached to 100 feet of 
log-line and sinker. This device measures the total drag 
induced on the line and sinker when it is towed from 
the stern of the launch. The instrument is calibrated 
to read in 1/10 knot units with a maximum range of 
zero to ten knots. The instrument is not accurate in 
the sense of conventional laboratory velocity meters. 
Accurate speed measurement at sea is extremely diffi- 
cult, particularly at low velocities. When suitably cali- 
brated, this instrument can be used to give a reliable 
estimate of ship speed. The instrument was calibrated 
periodically by determining the time of travel over a 
known distance at various constant launch speeds. 

The underway sampling equipment mounted on the 
launch ready for field use is shown in Figure 2. 


Rhodamine-B Techniques 


Feurestein and Selleck (2) recently published the 
results of an extensive series of tests designed to ex- 
amine critically the physical and chemical character- 
istics of several tracers commonly used in surface water 
investigations. They discuss the effects of physical 
adsorption, photo-chemical decay, temperature, pH, 
etc., on rhodamine-B, Pontacyl Pink-B, fluorescein and 
orzan. These authors concluded their report by stating 
that Pontacyl Pink-B is the most suitable of the tracers 
studied for use in natural surface waters. Rhodamine-B 
was found to be less suitable because of the unfavorable 
characteristic of adsorption on suspended solids indig- 
enous to San Francisco Bay. The magnitude of adsorp- 
tion was found to be large and highly variable with the 
chlorosity of the surface water, the higher the chlorosity 
the lower the adsorption. Both tracers were found to 
be significantly temperature sensitive, fluorescence de- 
creasing with increasing temperatures. pH values less 
than 5 or greater than 10 result in a decrease in fluores- 
cence for rhodamine-B. No such pH effect was demon- 
strated with Pontacyl Pink-B. Both tracers show slight 
photochemical decay rates up to about 30 hours ex- 
posure time, after which significant decay rates are 
evident for both, rhodamine-B decaying faster than 
Pontacyl Pink-B. Most of the above mentioned char- 
acteristics for rhodamine-B were discussed, at least in 
part, by Pritchard and Carpenter (8) in their study of 
diffusion in Chesapeake Bay, as outlined in a previous 
progress report. These authors, however, did not include 


a study on the adsorbtion characteristics of rhodamine-B 
in waters with high concentrations of suspended solids. 
The question that followed, then, was just how 
significant were the above mentioned characteristics for 
rhodamine-B when used in an oceanographic environ- 
ment as in this investigation? The pH of the surface 
waters near the Orange County outfall averaged from 
7.6 at the boil to 8.3 at non-boil stations. The raw 
rhodamine-B acetic acid mixture has a pH of approx- 
imately 3.5 and when mixed with methanol for specific 
gravity adjustment, the pH is approximately 4.0. The 
minimum pH observed from grab samples taken 30 
seconds after dye release was 4.2, this value increasing 
to 6.5 in grab samples taken 114 minutes later. There- 
fore, it appears that pH was no problem except possibly 
during the first few minutes following dye release. 

Temperature variations during field tests at Orange 
County were observed to be minimal. The maximum 
variation in surface water samples between boil and 
downstream stations was approximately 2°C. The tem- 
perature effect was also compensated for by periodic 
calibration of the fluorometers while following the dye 
patch. Before release of the dye, the instrument was 
zeroed by pumping local water through the unit. A 
standard solution was then prepared also using local 
water. With this standard (normally 1 ppm at the be- 
ginning of a test) the instrument sensitivity was set 
to give the desired dial reading, say 50 units, corre- 
sponding to the standard. This procedure was repeated 
at periodic intervals to check the initial standard and 
again when an increase in instrument sensitivity was 
required. 

Photochemical decay, assuming the values deter- 
mined by Feuerstein and Selleck apply to conditions 
existing in the ocean, did not appear to be of major 
concern in the studies conducted to date. The maximum 
daylight exposure time was seven hours. According to 
the above mentioned report, photochemical decay for 
rhodamine-B becomes serious after approximately 30 
hours of exposure. 

A more important factor to consider in studies around 
an existing outfall, such as Orange County, is the 
turbidity of the diluted waste found at or near the 
boil. Table 1 shows the results of a laboratory test 
using primary effluent from the city of Los Angeles’ 
Hyperion Treatment Plant tagged with rhodamine-B. 
The undiluted effluent shows a pronounced decrease 
in dial reading compared to the tap water blank. At a 
dilution (tap water to effluent) of 10:1, the effluent 
mixture reading equals that of the blank. This decrease 
in sensitivity is probably due to the light scattering 
characteristics of the particulate matter contained in 
the effluent. The initial dilution at the Orange County 
Outfall, as determined by salinity measurements, 
ranged from 7:1 to 20:1. Therefore, on the assumption 
that the Hyperion and Orange County Treatment 
Plant effluents are similar, the light scattering effect 
was assumed to be minimal. 

The problem of adsorbtion of rhodamine-B on sus- 
pended sediments similar to those found in San Fran- 
cisco Bay is probably not of major concern in the open 
ocean, the reason for this being that suspended sediment 
concentrations are generally negligible and chlorosity 
is high. A problem can exist, however, if studies are 
conducted near existing outfalls where the effluent 
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TABLE 1 


Hyperion Primary Effluent Fluorometer No. 2 

Suspended Solids = 90 ppm Flow Door Assembly 

Dilution with Tap Water Filters 

Temperature of Effluent Primary: 110-832 

and Water = 21°C Secondary: 110-833 
+110-823 


N.D. 1% 








Sample Dilution Dial reading Remarks 





Blank -- 64.5 1 ppm 
thodamine-B 





1 ppm 

rhodamine-B 

well mixed 
sample 


Effluent — 57.5 





1 ppm 

rhodamine-B 

well mixed 
sample 


Tap Water and 
Effluent 2:1 60.5 





1 ppm 

rhodamine-B 

well mixed 
sample 


Tap Water and 
Effluent 3:1 62.25 





1 ppm 
Tap Water and rhodamine-B 








Effluent 5:1 63.5 well mixed 
sample 
1 ppm 
Tap Water and rhodamine-B 
Effluent 10:1 64.5 well mixed 
sample 





All readings taken immediately after addition of dye. 


contains high concentrations of suspended sewage solids. 
Tables 2 and 3 show the results of adsorbtion studies 
conducted using Hyperion primary effluent. The values 
in Table 2 indicate the change in concentration with 
time for various dilutions of effluent, all containing 1 
ppm of rhodamine-B. There was an initial difference 
in fluorescence level for the undiluted effluent which 
was probably due to light scattering. In the other 
samples there was little if any initial difference. There 
appeared to be a negligible decrease in fluorescence for 
any of the samples for the eight hour test period. The 
values indicated in Table 3 are for a similar test using 
0.1 ppm rhodamine-B. There was again an immediate 


’ TABLE 2 


Fluorometer No. 2 
Non-Flow Door 
Cuvette No. 2 


Adsorbtion Study 

1 ppm rhodamine-B 
Hyperion Primary Effluent 
Suspended Solids = 90 ppm 
Dilution with Tap Water 


























Average 
Sample Dilution Time dial reading Remarks 
1103 52 6 et & 
Tap Water = 1118 51 
Blank 1400 51 
1923 52 
1107 50 Ti= 21.9°C 
Effluent — 1122 50 Floc on bottom 
1400 50.75 Heavy floc 
1923 50 Heavy floc 
1111 52 Ti= 22°C 
Tap Water 1123 62225 Floc on Bottom 
and Effluent 2s) 1400 51 Floc on Bottom 
1924 52.25 Floc on Bottom 
1116 51,5 T= 21.9°C 
Tap Water 1125 50.5 Fine floc on 
and Effluent 10:1 bottom 
1400 50.25 Fine floc on 
bottom 
1924 51.25 Fine floc on 
bottom 


All samples agitated for 10 minutes and allowed to settle before analysis. 
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initial decrease of fluorescence in all samples which ap- 
peared in this case greater in the samples of higher 
dilution. The important point is that only in the un- 
diluted effluent sample was there any appreciable 
decrease in fluorescence with time. 

The concentration of suspended solids in all water 
samples taken at the boil during this investigation 
ranged from 14.6 to 81.5 ppm. Because of this low 
range in solids concentration and the data obtained 
from the laboratory tests mentioned above, it appeared 
reasonable to assume that physical adsorbtion was not 
an appreciable factor in this investigation. In order to 
substantiate the results of the lab tests, a series of 
experiments were designed to compare rhodamine-B 
and Pontacyl Pink-B under actual field conditions. On 
May 28th and 29th, 1963, 5 lb. slugs of these tracers 
were released just downstream from the boil at the 
Orange County outfall site. Pontacyl Pink, used as a 
control, was released during the morning hours followed 
by rhodamine-B in the afternoon. Figure 3 shows the 
results obtained when the data from both days are 
superposed. The least-squares regression lines are also 
indicated with corresponding equations. Both equations 
fit the data very well, having correlation coefficients of 
0.98 and 0.96 for rhodamine-B and Pontacyl Pink 
respectively. The difference between the two curves 
can be attributed to many factors. If due to physical 
adsorbtion on suspended particles, it would have been 
necessary for a large quantity of rhodamine-B to have 
been removed initially, an occurrence that is not borne 
out by previous laboratory or field observations. The 
slopes of these curves appear to indicate that little, if 
any, adsorbtion took place with time. It is possible 
that the difference in these results can be attributed 
to varying oceanographic conditions existing in the 
vicinity of the boil at the time of release. Previous 
tests (as discussed later in this report) indicated that 
concentration vs time curves, for tests conducted in 
the same area during the same day, vary to an even 
larger degree than the curves shown in Figure 3. 

Another physical characteristic of both rhodamine-B 
and Pontacyl Pink-B that presented a problem in field 
sampling was the so-called concentration “quenching” 


TABLE 3 


Fluorometer No. 2 
Non-Flow Door 
Cuvette No. 2 


Adsorbtion Study 

100 ppb rhodamine-B 
Hyperion Primary Effluent 
Suspended Solids = 58.6 ppm 
Dilution with Tap Water 



































Average 
Sample Dilution Time dial reading Remarks 
0945 53.0 T= 21.3°C 
Tap Water — 1143 53.0 T = 21.38°C 
Blank 1712 52.5 T = 22°C 
0953 51.25 Pe 216 °C 
Effluent — 1134 51.0 T = 21.9°C 
1716 49.5 T = 22°C 
1720 50.75 re-mixed sample 
0959 49.75 T = 21.5°C 
Tap Water PA | 1139 49.75 T as 215°C 
and Effluent 1724 49.25 T = 22°C 
1735 49.25 re-mixed sample 
1003 49.5 T = 21.4°C 
Tap Water 10:1 1130 50.0 T = 21.4°C 
and Effluent 1730 49.5 T = 22°C 
1735 49.6 re-mixed sample 


All samples agitated for 10 minutes and allowed to settle before analysis. 
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effect. When either dye is present in high concentrations 
it tends to absorb the ultra violet exciting light used 
in the fluorometer. Therefore, dye concentration is not 
a linear function of the instrument dial reading. Figure 
4 is a typical calibration curve for rhodamine-B ob- 
tained when the instrument sensitivity was set so that 
1 ppm equals 50 dial units. The severe non-linear read- 
out characteristic is the result of the ‘‘quenching”’ effect. 

Analysis of strip chart data from continuous-flow 
measurements taken at the early phases of a test run 
required careful interpretation. The problem was to 
determine which portion of the calibration curve was 
to be used in data reduction. Figure 5 is an illustration 
of the strip chart data taken continuously along the 
longitudinal axis of a typical dye patch. As the con- 
centration increases the instrument senses the values 
indicated on the lower portion of the calibration curve 
in Figure 4. Near the center of the patch the concen- 
tration gradually increases to 1 ppm, and then beyond. 
The instrument now senses the upper region of the 
curve where increasing concentrations have lower 
apparent fluorescence levels; thus the readout indicates 
a drop which corresponds to the maximum concentra- 
tion. The subsequent decrease in concentration is 
indicated by an increase in dial reading back to approx- 
imately 1 ppm., after which the instrument again senses 
fluorescence levels corresponding to the lower portion 
of the calibration curve and decreasing concentrations 
are indicated by decreasing dial readings. In this ex- 
ample the maximum concentration is actually 3 ppm. 
instead of an apparent maximum of 1 ppm. 

It was observed that when using 5 and 10 Ib. slugs 
of rhodamine-B dye, it took from 10 to 40 minutes 
before the maximum concentration in the patch was 
lowered to less than 1 ppm. Therefore, during the early 
phases of a given test, grab samples together with 
continuous-flow samples were taken concurrently in 
order to establish which portion of the calibration curve 
was to be used in later analysis of the data. The grab 
samples were analyzed at the conclusion of the test by 
passing diluted aliquots through the fluorometer. 


3. DILUTION STUDIES 


Dilution studies were conducted in order to provide 
back-up information for phytoplankton productivity 
investigations that were conducted within the sewage 
field. The rate of dilution was estimated both by in- 
stantaneous release of slugs of dye and by continuous 
release from a skiff anchored near the surface boil at 
the Orange County outfall site. A summary of the 
oceanographic conditions, empirical dilution param- 
eters, and equations for the instantaneous release 
studies in both the waste field and in control areas 
was presented in reference (1) and is included here for 
further discussion (Tables 4 and 5). The results obtained 
from the continuous release experiments in the waste 
field are summarized in Tables 6 and 7. 

It should be pointed out that, with one exception, the 
results of both the instantaneous and continuous re- 
lease experiments, do not give the true rate of dilution 
of the entire waste field under investigation. The 
reason for this is the difference in size of the waste and 
dye fields. The rate of eddy diffusion has been shown 
to be a function of the scale of the turbulent eddies 


participating in the diffusion phenomenon. Therefore, 
unless the waste and dye fields in a test area are of 
comparable size, the rates of diffusion and therefore of 
dilution will not be equivalent. The exceptional case is 
the experiment conducted on 17 June 1964 when the 
entire waste field was tagged by continuous injection of 
dye at the Orange County treatment plant on shore. 

Regardless of these difficulties, the data presented do 
give a measure or indication of the rate of dilution 
occurring on any given day. The results of individual 
experiments shown in Tables 5 and 7 can be compared 
to determine relative rates of dilution. In this respect 
these data can be used to estimate and compare the 
influence of physical dilution on primary phyto- 
plankton productivity within the waste field. 


4. MATHEMATICAL MODELS—RELATIVE 
DIFFUSION 


The information discussed in this section was given 
partial treatment in the 1962-63 Annual Report (1) 
and was also included in a paper (4) presented at the 
1964 meeting at Monaco of the Commission Interna- 
tionale pour |’Exploration Scientifique de la Mer 
Mediterranee (Tibby, Foxworthy, Oguri, and Fay, 
“The Diffusion of Wastes in Open Coastal Waters and 
their Effects on Primary Biological Productivity.”’) 
Some of the material is presented again at this time to 
serve as a basis for the discussion and conclusions which 
follow. 

During the first few field experiments it was observed 
that the dye patches elongated in the direction of the 
surface current. This suggested that, neglecting the 
possible effect of a vertical velocity gradient near the 
surface, the process was anisotropic, i.e., diffusing at 
different rates in the horizontal plane. The highly posi- 
tive density gradient normally found within the surface 
waste field from the Orange County outfall also indi- 
cated that the rate of vertical diffusion was much less 
than the horizontal. Under these conditions the diffusion 
process was considered to be of a three-dimensional 
anisotropic type. This assumption will be further dis- 
cussed in a later section. 

The three-dimensional statistical model proposed by 
Gifford (5) and extended by Okubo (6) to cover the 
case of relative diffusion in a homogeneous and station- 
ary turbulent field was considered to be directly ap- 
plicable to the case in point. This model, as discussed 
in reference (1), is based on the assumptions that 
diffusion proceeds (independently) at different rates in 
the horizontal and vertical planes and that the spatial 
distribution of diffusing substance is essentially Gaus- 
sian in all directions. The later assumption, assumed 
to be valid for the average concentration distribution 
of a large number of dye releases diffusing under similar 
oceanographic conditions, leads to the following equa- 
tion: 


C(x,y,z,t) (1) 


M x? y? z2 
 wa/ Der [a x22? |! /? ea =f 2a? a4 Qzey? T al 





where C is the average concentration at a point 2, y, z 
and time t, M is the amount of dye initially discharged 
from an instantaneous point source, and o? are the 
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average values of the variances of the concentration 
distribution in the coordinate directions as a function 
of time. Gifford also points out that no a priorz assump- 
tion about the form of the variances is necessary to 
make use of equation (1), the only stipulations being 
those given above. 

An inspection of equation (1) indicates that the 
maximum concentration as a function of time is given 
by the first term on the right, i.e., by setting x = y = 
z = O equation (1) becomes: 


M 
wr/ 2m [a2 25,2 }! /2 (2) 

Equation (2) is of special interest because it is appli- 
cable to the disposal of a patch of waste material in the 
sea where the decrease in maximum concentration with 
time is of primary interest. It is relatively easy to obtain 
field data for the verification of equation (2), the prob- 
lem centering on determining the variances and maxi- 
mum concentration within the patch as a function of 
dispersion time. 

As already stated, the dye patches elongated with 
time in the direction of the surface current. In most 
experiments this elongation was quite marked with the 
overall shape resembling an ellipse or, considering the 
three-dimensional picture, an ellipsoid. 

Equation (1) can be shown to define a series of three- 
dimensional ellipsoidal surfaces of equal concentration, 
the major and minor semi-axes being defined (Dia- 
chischin (7) ) by the following equations: 


nee = 











a? = In (SE) " 3 (3) 
bt = In (S")' (4) 
ln (Se) 32 (5) 


where a, 6, and c are the semi-axes at any given time 
defining a surface with equal concentration C’, and 
Cmax is the maximum concentration at the center of 
dye mass. By repeated underway sampling along the 
principal axes of the patch, it was possible to determine 
the distribution of the corresponding dye concentra- 
tions. From these data the average values of a, }, c, 
C’ and Cmax were determined, and the values of the 
coordinate variances then estimated by use of the above 
equations. 

Figures 6 and 7 show a series of vertical photographs 
taken during the initial phases of an instantaneous 
release experiment. Five pounds of rhodamine-B were 
released downstream from the boil (lower center of 
Figure 6). The time-variation of the size and shape of 
the patch can easily be seen. The patch elongated in 
the direction of motion, and the extent to which the 
shape of the patch approached that of a true ellipse 
should be noted. 


5. SAMPLING PROCEDURES 
The sampling procedure followed in all dye patch 
studies is summarized as follows: 
1. The fluorometer was calibrated with standard 
solutions of tracer dye. The initial sensitivity setting 


was 1 ppm equal to 50 or 75 fluorometer units. The 
instrument was “zeroed” using water from the waste 
field as a blank. 

2. A dye tag was initiated by release of 5 or 10 
pounds of dye solution with specific gravity of approxi- 
mately 1.02. The release was usually completed in less 
than one minute. 

3. During the first few minutes following dye release, 
grab samples were taken at the point of apparent max- 
imum concentration. 

4. Continuous flow sampling commenced in the 
center of the patch to establish the depth to maximum 
concentration. 

5. With the underway sampler set at this depth, 
continuous flow samples were taken periodically along 
the principal horizontal axis of the patch. Dye con- 
centrations were recorded continuously on a strip chart 
readout attached to the fluorometer. 

6. Continuous flow depth samples were taken peri- 
odically from the research launch while located in the 
zone of (visual) maximum concentration. 

7. The fluorometer was recalibrated at approximately 
one hour intervals. 

8. The vertical distribution of temperature and 
salinity was determined at various time periods by 
personnel on the R.V./Velero IV. 

9. Aerial photographs were taken at regular time 
intervals to aid in subsequent data reduction and 
analysis. 


6. EXPERIMENTAL RESULTS 


Figures 8 to 10 are examples of the measured values 
of the concentration distribution along the principal 
x, y and z-axes. These data are for single crossings of 
the patch at various indicated average times after dye 
release. The observed concentration distributions have 
been fitted with hypothetical Gaussian curves obtained 
from the estimated values of the variances and meas- 
ured values of Cmax. It can be seen that the computed 
normal curves fit the observed data reasonably well 
during the early phases of diffusion. At later times the 
fit is not nearly as close. As described above, better 
results could be obtained if averages of several crossings 
taken at the same time could be used to establish the 
observed points and to estimate the magnitudes of the 
variances. Unfortunately, data of this type, at least for 
studies of relative diffusion, are almost impossible to 
obtain using the present sampling techniques. However, 
the values obtained for the variances by this method 
are believed to be a reasonably good estimate of the 
average values. 

It should be mentioned that the spatial concentration 
distributions shown in Figures 8 to 10 are neither the 
best nor the poorest recorded. In some experiments the 
distributions were skewed to the right or left of Cmax. 
In these cases the basic assumption of a Gaussian dis- 
tribution breaks down. However, on the basis of the 
majority of test results, it appears reasonable to con- 
sider a normal distribution as a useful, if not completely 
accurate working rule. 

Once the values of the diffusion parameters were 
determined, it was possible to ascertain the empirical 
relationships between these parameters and diffusion 
time. It was found that when the variances were plotted 
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against time on log-log paper, the points in most cases 
defined a straight line. Empirical relationships for the 
variances of the form: 

o? == k(T)® (6) 
were then computed by the method of least-squares. 

Figures 11 through 13 indicate the results obtained 
from the study conducted at the Orange County site 
on 11 April 1963. For this experiment five pounds of 
dye solution was discharged at a point about 50 yards 
downstream from the surface boil. The dye patch was 
sampled at frequent time intervals for a period of about 
6 hours. 

The data as shown in Figures 11, 12 and 13 have 
been fitted with the appropriate least-squares trend 
lines and the empirical equations for the variances are 
indicated. Visual inspection of these figures shows that 
the correlation between o*,,, and T is quite strong. 
The o?, vs. T plot shows considerable scatter about the 
trend line, indicating that the correlation is less sub- 
stantial. Data from other patch experiments give 
similar results, although the o?, vs. T trend line in most 
cases is much better than shown here. 

The most direct method of testing the validity of 
equation (1) is to work directly with the observed 
values of the maximum dye concentration within the 
patch, that is, by use of equation (2). Substituting the 
empirical expressions for the variances into the equation 
for Cmax yields: 


Cimax = 0.037 (T)-1-*! (Ib/ft?) 


Taking the specific weight of sea water as 64 lb/ft, we 
can also express the results in units of parts of dye to 
million parts of sea water (ppm) as: 

0.037 


= ——— 6 SUN 
Cmax 64 x 106(T) 


= 578(T)~'-*!(ppm) 


This expression can be compared to the measured 
values of Cmax as a function of dispersion time taken 
from the fluorometer record. The test of the validity of 
equation (2), the field sampling procedure, and of the 
analytical method used to compute the variances, is 
the agreement between these two expressions for Cmax. 

The comparison between the measured and com- 
puted values of Cx is shown graphically in Figure 14. 
The circled points are values of Cnax taken from the 
fluorometer record. These points have been fitted with 
a least-squares trend line as shown. The expression for 
Cmax derived by equation (2) is indicated by the dotted 
line. The two curves, for this particular set of data, 
agree quite favorably. 

It is important to note that the values of Cnax used 
in determining the variances by use of equations (3) 
through (5) do not necessarily correspond to the values 
of Cmax plotted as ‘field data.’’ The individual maximum 
concentration for a particular series of transects along 
the z, y, and z axes (as shown in Figures 8 through 10) 
were used to compute o*, y,,. The value of Cmax selected 
for the ‘field data” plot corresponds to the highest 
value of Cnax obtained from this, or even a later, series 
of samples. 

Attention should also be called to the two-dimensional 
form of equation (2) in which the vertical component 


of diffusion is neglected and diffusion is assumed to 
take place at different rates in the horizontal plane. 
Assuming the dye is concentrated in a surface layer of 
d thickness, the two-dimensional form of equation (2) 
is 

M/d 


or [Fx2ay? }tv2 


Coe = (2a) 


where d is the average depth of the dye layer. Using 
the estimated values of 6? and ¢,? as before, and an 
average depth of the dye layer of 10 feet, we obtain: 


aie: 5/10 
max Qn [68.8 X 13.5]? 


= 0.00275(T)-!-55(1b/ft) 
= 43(T)~'>(ppm) 


The above equation is also plotted in Figure 14 for 
purposes of comparison. The effect of neglecting the 
vertical component of diffusion is obvious. The values 
of Cmax Obtained using the two-dimensional model, for 
these data, do not furnish nearly as accurate estimates 
as the three-dimensional one. 

Several mathematical models recently have been 
proposed for the case of horizontal diffusion from an 
instantaneous vertical line source in an isotropic, hom- 
ogeneous and stationary turbulence field in the sea. 
These models indicate that Cnax is proportional to 
(T)-, where the exponent n takes on values of 2 or 3. 
The data summarized in Table 5 indicate that in these 
experiments the value of the exponent in the Cmax 
equations ranged from 1.53 to 2.29. In no case did n 
approach a value of 3. This suggests that the horizontal 
diffusion models having the form Cmax « (T)—’, as pro- 
posed by Joseph and Sender (8), Okubo and Pritchard 
(9), and Schonfeld (10), might be applicable to some 
of the present experimental data. The model proposed 
by Joseph and Sender has been examined in a previous 
Progress Report (11) and will be used here as a basis of 
discussion. Their equation which is reportedly valid for 
eddy scales of 10 to 1500 kilometers, takes the form 


Cien) = pes ex ~ (Ge) 8 


where m’ is the amount of dye discharged per foot of 
depth, P is the ‘‘most probable velocity of diffusion” 
assumed equal in all directions in the horizontal plane, 
and r is the radius from the center of dye mass to a 
point having a concentration C. 

Equation (8) can be solved in many ways (Okubo 
(6) ). The method proposed by Ichiye (12) appeared to 
be the most manageable and was applied to several 
sets of data obtained during the first series of patch 
experiments. The results obtained for the data collected 
on 17 August 1962 are typical. The observed dimensions 
of the 0.05 ppm contour were plotted at various time 
periods, and the corresponding horizontal areas were 
measured with a planimeter. The radius r of an ‘‘equiv- 
alent circular area’? was then determined. Using esti- 
mated values of the ‘‘diffusion velocity” P in equation 
(8), it was possible to compute theoretical values of r 
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for different diffusion times and thus to fit a theoretical 
curve to the observed points. Figure 15 shows the re- 
sults obtained by this method of analysis. A diffusion 
velocity of 1.7 ft/min. appeared to give the best com- 
parison between observed and theoretical curves. The 
value of P for this test agrees favorably with those 
determined for other experiments, as well as with those 
given by Joseph and Sender in their original paper (8) 
for diffusion on a much larger scale. 

Assuming a value of P equal to that determined 
above, it is possible to compare the observed and theo- 
retical values of Cnax for the data of 11 April. Assuming 
an average surface dye layer thickness of 10 ft. as be- 
fore, and setting r = 0 in equation (8), we have: 


m’ 


Gos = 2 2PHt Oi 


429(T)-2 


The results of this analysis are also shown in Figure 14. 
It can be seen that for this particular experiment the 
model proposed by Joseph and Sender gives results that 
compare reasonably well with the observed maximum 
concentrations. It must be emphasized that models 
similar to equation (8) and (9) are derived on the as- 
sumption that the turbulence is horizontally isotropic. 
Therefore, these models cannot be expécted to describe 
accurately the spatial distribution of dye in which the 
turbulence is wholly or in part anisotropic. This can 
best be seen in Figure 16 where the actual and estimated 
dye distribution is shown for the data of 17 August 
approximately 110 minutes after dye release. 

In connection with the models just described, it is 
interesting to speculate further on the assumption of 
horizontal anisotropy. Ichye (13) reports that the 
elongation of dye patches can be the result of both 
anisotropic conditions as well as “shear effect’? due to 
the existence of a velocity gradient with depth. He 
states that photographs of dye patches were observed 
to form a narrow ‘“‘tail’’ of dye, at depth, which he 
attributes to a vertical velocity differential. In the ex- 
periments conducted by the Allan Hancock Foundation, 
the underway sampler normally was set at a depth of 
0 to a maximum of about 2 feet. At those shallow depths 
(which can be considered essentially surface sampling) 
the length of the dye patch was usually much greater 
than the width. This is shown in Figures 8 and 9. This 
fact does not preclude the effects of a vertical velocity 
gradient in the thin upper layers of the water surface, 
but at least gives some substance to the assumption of 
horizontal anisotropy. If shear is important, its effect 
is included in the measured values of o7,,, shown in 
Table 5. Field measurements and analytical techniques 
devised to separate the shear effect from true anisotropy 
must be left for future study. 


7. MASS BALANCE ANALYSIS 


In any study of relative diffusion it is necessary 
periodically to re-estimate the amount of dye remaining 
within a patch. An analysis of this type is termed a 
“mass or materials balance.’”” Depending on the accu- 
racy of the calculations, a mass balance may indicate 
the amount of dye lost through physical adsorption, 


photochemical decay, or other factors. In Pearson’s (14) 
review of tracer methodology, he points out the im- 
portance of such estimates in diffusion studies utilizing 
artificial tracers. 

In the present experiments the amount of dye was re- 
estimated by calculating M in equation (2) using ob- 
served values of Cmax and the empirical values of the 
variances at different time periods. Thus we have: 


M’ = CpaxtV/2r [525 y25,? ]/ (10) 


Where M (referred to as M’) is in pounds and the units 
of Cmax are lb/ft. Obviously M’ computed by this 
method is a function of the diffusion time as long as 
the slopes of the curves for the observed and calculated 
values of Cmax are not parallel. Therefore, for purposes 
of comparison, the average value of the diffusion time 
was used in all mass balance calculations. Applying 
this method of analysis to the data collected on 11 
April yields: 


M’ = [4.66 X 10-4] (684(T)—!-® ] [17.2(T)?-9!] 
OL Os 

At T = 150 minutes: 
M’ = 6.08 (Ib) = 122% recovery 


It is seen that for this experiment M’ > M, or that 
more dye was recovered at 150 minutes than was orig- 
inally released. The discrepancy possibly can be ac- 
counted for by induced errors in sampling technique 
or by increased indicated fluorescence levels due to the 
presence of particulate matter indigenous to the waters 
within a waste field. 


Table 5 summarizes the calculated values of M’ for 
all dye patch studies. The indicated values range from 
40 to 244 percent. Bowles et al. (15), in their studies 
of the diffusion characteristics of fluorescein in the 
English Channel, reported mass balances ranging from 
35 to 304 percent. They speculated that the higher 
values of M’ were probably due to the presence of 
particulate matter. Pritchard and Carpenter (8), in 
their study of diffusion in Chesapeake Bay, reported 
that mass balances of 80 to 90 percent were realized 
using rhodamine B as a tracer. 

It is interesting to note that the higher values of M’ 
indicated in Table 5 correspond to experiments carried 
out within the waste field. The higher concentration of 
particulate matter in these waters could have been re- 
sponsible for the higher indicated recoveries. However, 
when one considers the initial dilution at the boil (values 
for these experiments ranged from 7/1 to 20/1.) plus 
the rapid dilution as the field moves downstream, it is 
not likely that particulate matter was a serious prob- 
lem. This is confirmed by transparency measurements 
taken in the field. 


The values of M’ for experiments conducted in the 
background sites outside the waste field ranged from 37 
to 78 percent. These waters were characterized by rela- 
tively low natural turbidity. In all but one of these 
experiments the stabilities were much lower than that 
observed in the waste field at Orange County. The rate 
of vertical diffusion, with one exception, was observed 
to be much greater in these areas. Deeper dye penetra- 
tion made vertical sampling more difficult, and it is 
possible that the observed data were less accurate for 


26 


OUR 
S) 
° 


RADIUS OF CONT 


rr re) 
(eo) 


o 


FATE OF WASTES IN THE MARINE ENVIRONMENT 


RADIUS VS. TIME 
FOR 
0.05 PPM CONTOUR 
ORANGE COUNTY 
I7 AUGUST 1962 


OBSERVED DATA 


THEORETICAL CURVE 
J&S, P= 1.7 Ft./Min. 


20 40 60 80 100 120 140 
TIME - (MINUTES) 





FIGURE 15 


FATE OF WASTES IN THE MARINE ENVIRONMENT 27 


0.05 PPM CONTOUR 
AT 

T = 1O MINUTES 

ORANGE COUNTY 

I7 AUGUST 1962 


0.05 PPM CONTOUR 


ree 
re 
Oo 
<9) 
o 


EQUIVALENT AREA 


Fe 
- 
= 


r= 146 FT. 





FIGURE 16 


28 FATE OF WASTES IN THE MARINE ENVIRONMENT 


these experiments. This possibility is substantiated 
somewhat by considering the data collected on 15 
August 1963 near the Redondo Canyon in Santa Monica 
Bay. On this day the water was particularly clear, but 
the average stability was comparable to values obtained 
within the waste field at Orange County. Vertical diffu- 
sion was suppressed as indicated by the empirical 
equation for the diffusion parameter o?, (Table 5). 
However, the value of M’ (78 percent) compares favor- 
ably with those determined within the waste field. It 
appears that on days of relatively low stability which 
results in an increase in vertical diffusion (see next 
section), the sampling technique may be less accurate. 

Whatever the reason for the discrepancy in the cal- 
culated values of M’ within and outside the waste field, 
it appears that it was not the result of photochemical 
decay, water temperature fluctuations, or physical 
adsorption. All experiments were conducted at or very 
near the water surface, where photochemical and tem- 
perature effects presumably would be similar regardless 
of the location of the test area. Physical adsorption, if 
appreciable, probably would occur more readily within 
the waste field where the concentrations of organic 
solids were higher than at the background sites. If any- 
thing, lower values of M’ would be expected for ex- 
periments conducted within the waste field but just 
the reverse was true in these experiments. 


8. EFFECTS OF WIND AND STABILITY 
ON DIFFUSION PARAMETERS 


One of the principal objectives of this research has 
been to demonstrate, if possible, the effect of water 
column stability and wind speed on the rate of vertical 
eddy diffusion. For the case of relative diffusion at the 
surface of the ocean, it has often been speculated (16) 
that the presence of a strong positive density gradient 
with depth (a highly stable water column) can only 
result in a supression of vertical turbulence and a con- 
comittant decrease in the rate of vertical eddy diffusion. 
This speculation has been confirmed somewhat by 
Kellogg (17) whose investigations on the diffusion of 
smoke puffs in the atmosphere showed that a decrease 
in the mass rate of puff growth occurred with increasing 
atmospheric stability. On the basis of Kellogg’s work 
it was believed a similar effect could be demonstrated 
in the ocean. 

To the oceanographer the term “stability” generally 
refers to the definition presented by Sverdrup et al. (16): 
Op es tool tas0P 


Ee ee AZ ned (11) 


where p is the density of sea water defined by ps,6,p 
that is, p is the density zn sctu at a salinity s, a tempera- 
ture @ and a pressure p. A distinction is usually made 
between the density in situ and that at atmospheric 
pressure, ps,e,o. The difference between the two is due 
to the compressibility of water at the particular tem- 
perature and salinity. When working in relatively 
shallow depths the difference may be neglected. 

In order to avoid writing a large number of decimals, 
the density (ps,6,0) is normally expressed by the symbol 
sigma-t, defined by: 


Ca (ps,6,0 = 1) 10? (12) 


An approximate form of (11), accurate within the first 
300 feet of the water surface, is expressed as: 


, _ dor ag 
= dz Pat (13) 


Where Z is the depth in meters. Actually E’ is the slope 
of the o; curve in relation to depth and can be positive, 
neutral, or negative. A water column is classified as 
stable, indifferent, or unstable depending on the sign 
of EK’. In these experiments stability has been arbi- 
trarily defined, for reasons of convenience, as 

S’ = EB’ Xx 10? 

The importance of selecting the proper depth interval 
when computing and interpreting values of S’ within 
a waste field has been discussed by Gunnerson (18). 
Early experiments conducted around the Orange County 
outfall indicated that at 50 to 100 feet from the axis of 
the boil the depth of the field was less than 3 meters. 
Consequently, temperature and salinity data for sta- 
bility computations were collected at 0, 3, 6 and 12 
meter depths. Stability data for the 0-3 meter interval 
were found to give the strongest correlation with mea- 
surements of vertical diffusion and all values of S’ used 
in this report therefore are for this depth interval. 

Stability within the waste field was highly variable. 
S’ ranged from highly positive values of 500 to 5800 
near the boil, to low positive and sometimes negative 
values at 200 to 250 minutes after dye release. Figure 
17 is a semi-log plot of the observed values of 8’ de- 
termined on 11 April 1963. It is evident that the sta- 
bility decreases quite rapidly from an initial high 
positive value at the boil to values of S’ approaching 
that of the background water at around 300 minutes 
downstream. Data collected during two experiments 
off Catalina Island are also shown for comparison in 
Figure 17. On these days the stability also varied with 
time, probably due to heating of the surface layers, but 
the individual values of S’ were considerably lower than 
corresponding values observed within the waste field. 

The space-time variation of 8’ within the field is the 
result of further dilution of the effluent mixture with 
adjacent and underlying water masses. As the field di- 
lutes with time, the surface salinity increases, thus 
decreasing observed values of S’. The same space-time 
variation has been observed with respect to inorganic 
nutrient and tracer dye concentrations as shown in 
Figure 18. Because of this pronounced time variation 
of stability within the field, it is necessary to use a 
time-averaging process to arrive at comparable values 
of S’. The average stability at any time is defined as 


ae 


Yi) .-194 (2S) aT + tie] 
wad aa Sieg 2 


where S’ is the average stability at time T minutes 
after dye discharge. 

For experiments where no boil samples were taken, 
the average of all observed boil stabilities (4500) was 
arbitrarily used. The computed average values of S’ 
for all experiments at or about 200 minutes after dye 
release are indicated in Tables 4 and 6. 
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Tibby and Barnard (19) pointed out that for a sur- 
face waste field, the magnitude of E’ is proportional to 
the amount of work required to mix the surface and 
underlying water masses. In Southern California, the 
energy required to promote vertical mixing in a stable 
surface layer usually is derived from wind stresses at 
the water surface. Another contributing factor related 
to wind action, is the breaking of wind-generated waves 
as reported by Johnson (20). 

From this discussion one would expect that the rate 
of dilution of a dye patch, as brought about by the 
horizontal and vertical diffusion processes, should be 
related to the average stability of the water column 
and average wind speed near the surface. 

In order to proceed with the investigation it was 
desired to develop a form of a dilution equation for 
the dye patch experiments. The ratio M/Cmax has di- 
mensions of 1°. Therefore the minimum patch dilution 
as a function of time was defined (Allan Hancock 
Foundation (1) ) as: 


M 
aha 





Dain = (ft) (14) 


Where Cmax is in Ib/ft* units taken from the empirical 
relationship for the observed field data. The dilution 
can also be expressed in terms of the three-dimensional 
diffusion parameters by substituting the value of Cmax 
given by equation (2), into (14) above 


Dan — w/ lr [F265 252" jae (15) 
The dilution equation for the data collected on 11 April 
1963 is: 


Dain = = 114(T)!-*(ft?) 


5) 
0.438(T)—1-89 


The empirical dilution equations for all dye patch 
studies that were discussed earlier are summarized in 
Table 5. 

In order to show the effect of wind speed and water 
column stability on the rate of dilution, graphs of 
dilution vs. average wind speed were constructed at 
arbitrarily selected values of diffusion time. Figures 19 
and 20 show the results obtained at diffusion times of 
50 and 200 minutes within the waste field. Values of 
the average stability at or about 200 minutes (Figure 
20) are indicated above each point and the correspond- 
ing test numbers (Table 4) are indicated below each 
point, all figures beginning with Figure 18. Inspection 
of these Figures indicates that values of dilution tend 
to increase with increasing average wind speed. The 
slope of the visual trend line becomes steeper with 
time. The effect of stability on the dilution-wind speed 
relationship is partially indicated in Figure 20, where 
the plotted points appear to form two groups. The 
lower group of points contain the higher values of 8’. 
At any given average value of wind speed, lower values 
of dilution are associated with higher values of average 
stability. Inspection of equation (15) indicates that 
this effect is controlled by the relationship between S’ 
and the vertical diffusion parameter ¢,”, provided the 
horizontal parameters are unaffected by stability. 
Sverdrup et al. (16) not only speculated on the supres- 


sion of vertical turbulence by increasing stability, but 
also stated that horizontal mixing would, for the same 
conditions, be relatively unaffected, the reason being 
that horizontal mixing supposedly takes place along 
planes of essentially equal density which are, of course, 
not altered by changes in the vertical density structure. 

The influence of wind speed and stability on the 
individual diffusion parameters can be shown by the 
method of analysis used above. Values of the coordinate 
variances at 50 minute intervals of diffusion time were 
computed using the empirical equations shown in 
Table 5. The variances were then plotted against. the 
corresponding average wind speeds with the average 
value of stability indicated for each point. Figure 21 
through 24 are semi-log plots illustrating the results of 
this analysis 200 minutes after dye release. 

Figure 21 shows that, for the data collected within 
the waste field, there is a strong tendency for @,? to 
increase with increasing values of W. The scatter in the 
individual observations can, in part, be attributed to 
variations in stability. These data also appear to show 
that a sudden increase in values of @,? takes place when 
the average wind speed exceeds about 8 knots. This 
effect is even more pronounced on similar plots for the 
earlier time intervals (not shown). The apparent sudden 
increase in slope of the «2 — W relationship was dis- 
cussed by the Allan Hancock Foundation (21) where 
it was hypothesized that a “critical wind speed” might 
exist between 6 to 7 knots. The critical wind speed 
has been defined by Johnson (20) as the minimum wind 
speed at which wind waves break and at which vertical 
turbulence is greatly increased. Johnson summarized 
all previous work on this subject, and reported that 
existing evidence indicates that a critical speed is 
reached at about force 4 on the Beaufort Scale (11 to 
16 knots). 

The wind data for these experiments was obtained 
from two recording anemometer stations. One ane- 
mometer is mounted on the R.V./VELERO IV about 
30 feet above the water surface. Since the research 
vessel was underway at variable speeds and directions 
most of the time during these experiments, it was nec- 
essary to rely mainly on data collected at the Hunting- 
ton Beach Steam Station of the Southern California 
Edison Company, where a fixed recording type ane- 
mometer is mounted about 25 feet above ground on a 
pole placed near the beach. The instrument is located 
about 114 miles northwest of the terminus of the Orange 
County outfall. 

Because of variation in drag force, the wind speed 
at a given altitude over land is more variable and 
generally less than the corresponding speed over water. 
The exact relationship is complicated and it is usually 
necessary to resort to empirical methods to establish 
the over-water/over-land speed ratios. A recent study 
by Saville, McClendon and Cochran (22) indicates that 
the ratio of overwater to overland velocities at a height 
of 25-30 feet, for a fetch of 8 miles is about 1.31. This 
value was determined by simultaneous measurements 
of overland-overwater wind velocities near a large in- 
land water supply reservoir. An attempt was made to 
establish a similar relationship for the Orange County 
test area. Airspeed recordings for the same period from 
both anemometer stations were compared. These data, 
collected on 14 different days, cover the period from 
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August to December 1962. For this purpose, the record 
from the anemometer located on the Velero IV was 
taken when the vessel was either anchored or drifting 
slowly in the test area. The average overwater to over- 
land ratio was found to be 1.3 with a median value of 
1.18, also at a height of 25-30 feet. The results obviously 
compare favorably with those mentioned above. The 
fetch distances for these data were estimated to be at 
least 20 miles, with the wind moving in from the sea. 


Using an average ratio of 1.3, the apparent critical 
wind speed in the Orange County test area is estimated 
to be around 9 to 12 knots. This value compares favor- 
ably with that given by Johnson. The foregoing dis- 
cussion must be classified as interesting speculation 
until more data can be obtained, especially at higher 
wind speeds. 


The effect of stability on ¢,? is shown in Figures 21 
and 22 in which the background stations well outside 
the waste field have been added for comparison. Tak- 
ing, for example, an average wind speed of 8-10 knots 
it can be seen that for the sewage field data, points with 
higher average S’ show lower values of ¢,?. There are 
some exceptions, but in general the same tendency is 
indicated for other wind speeds. The background sta- 
tions, with one exception, show much lower average 
stabilities and correspondingly higher o,? than those 
within the field. This is to be expected in view of the 
previous discussion. The one exception is found in 
Figure 22 where the background station ‘‘e’”’ represents 
information collected over the Redondo Canyon on 15 
August 1963. The stability conditions on that day were 
relatively high, giving an average 8’ at 200 minutes 
very similar to high stability days within the waste 
field. The data for this particular background experi- 
ment fits well with those taken under similar conditions 
within the waste field. 


Some caution must be observed when comparing 
data from the Orange County stations during the initial 
period following discharge, and the background sta- 
tions. The area at and near the boil is characterized by 
extremely turbulent conditions. Also, the kinetic and 
potential energy of the effluent jet as it reaches the 
surface are partially converted to hydraulic head and 
the resulting pressure gradient promotes a radial spread 
of the waste in this area (see later Figures 26, 27 and 
28). It is not known just how far this effect is manifest 
on the downstream side of the boil but observation 
indicates that the distance may be as much as 1000 
feet. The horizontal and, to a lesser extent, the vertical 
dispersion of the dye may also be influenced by the 
presence of this induced pressure gradient. The influ- 
ence may have been appreciable even though the scale 
of the dye patch was much less than the simultaneous 
scale of the waste field. 


The relationship between g,”, stability, and wind 
speed can also be shown by plotting the computed 
values of the variance at different times against average 
values of S’. Values of @,? for wind speeds ranging from 
8 to 11.5 knots and from 0 to 8 knots are plotted against 
the average stability in Figure 25. Similar data from 
the background stations are also shown for comparison. 
The Figure indicates that within both ranges of wind 
speed, there is a tendency for @,? to decrease with in- 
creasing stability. The plot also indicates that for any 


given stability, higher values of ¢,? occur at average 
winds of 8-11.5 knots than at lower wind velocities. 


The number of measurements are too few to permit 
firm conclusions at this time but a trend toward a 
decrease in vertical diffusion with increasing stability 
does appear to exist as predicted by Sverdrup et al. 
The fact that these data also indicate consistently 
higher values of &,? at a given 8’, for average wind 
speeds between 8-11.5 knots than between 0-8 knots 
also adds to our confidence in the idea that there is a 
critical wind speed above which the rate of vertical 
diffusion increases sharply. 


A similar method of analysis was employed in at- 
tempting to establish a comparable relationship be- 
tween the horizontal variances, wind speed and sta- 
bility. Figures 23 and 24 are semi-log plots indicating 
the results of the analysis at a diffusion time of 200 
minutes. The plot of cy? against W for the waste field 
studies, indicates considerable scatter between indi- 
vidual points. There appears to be a trend toward 
decreasing cy? values with increasing wind speed, al- 
though there are insufficient data to warrant any solid 
conclusions along this line. Stability appears to have 
little or no definite effect on cy?. Although data from 
background areas do indicate that much lower values 
of cy” are associated with lower values of S’, the care 
which must be exercised in comparing data obtained 
inside and outside of the waste field has already been 
discussed. This is especially true of the lateral and 
longitudinal components of dispersion. 


Figure 24 shows the relationship between o,? and W. 
Again there is considerable scatter of the plotted points, 
but the data from within the field do indicate a trend 
of increasing o,” with increasing wind speed. As in the 
case of lateral dispersion, there appears to be no direct 
correlation between stability and o,2. Taking W between 
8 and 11 knots in Figure 24, it is seen that the point 
representing the highest value of S’ (846) does not lie 
at the highest value of ,? but is intermediate in position 
with respect to stability. In the range of 2 to 8 knots 
an increase in ox”? with decreasing 8’ is suggested but 
there are insufficient data on which to base a firm 
conclusion. 


9. THEORETICAL CONSIDERATIONS . 


Gifford (23) reviewed the theoretical basis for relative 
diffusion of a cloud of particles in atmospheric turbu- 
lence, including the earlier work of Brier (24) and 
Batchelor (27, 29) in which, on the basis of the similarity 
theory of homogeneous turbulence, certain dimensional 
predictions were proposed regarding the mean-square 
dispersion of two fluid particles. These predictions were 
reportedly valid for the so-called inertial subrange of 
turbulent eddy sizes where eddies are small compared 
to the large scale energy-bearing eddies of the mean 
motion. The existence of an inertial subrange implies 
that local isotropic conditions prevail, and this in turn 
depends on the state of the turbulence. If the Reynolds 
number of the turbulence is large enough, the turbu- 
lence should be “fully-developed’”’ and an inertial sub- 
range will exist. Okubo (9) reported that in most oceanic 
areas the Reynolds number of turbulence is sufficiently 
high to meet this requirement. 
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Batchelor (25) pointed out that when the distance 
between two diffusing particles lies in the inertial sub- 
range, the rate of change of the particle-separation 
depends on the initial separation lo, the diffusion time, 
and the parameters of the turbulence (the eddy vis- 
cosity and the rate of energy dissipation per unit mass, 
e.). If the lower limit of the inertial subrange in which 
the effects of viscosity are felt is small, Batchelor’s 
dimensional predictions yield: 


dé 2/3 cc ” 
qe (eo)?/*t, (for “small” t) (16) 


where ¢ is the mean-square dispersion at time ¢ in a 

eae direction (the projection of J on the 
£ axis). 

When the diffusion time is large enough, the relative 

motion of the particles no longer depends on &, and 
dé? 


ae o et?, (ti K t) (17) 


where t, is the time beyond which the initial separation 


2 


has no affect on « and is given by: 


ty = 1p2/8 el/8 (18) 


Integration of (16) and (17) for a constant rate of 
energy dissipation yields 
E(t) — &(0) « t? (16a) 
f(t) « t3 (17a) 


It is important to note, as Batchelor (25) further 
indicates, that as the diffusion time becomes very long 
the two particles wander independently, and the com- 
ponents of the relative dispersion of the two particles 
reduce to twice that of a single particle. In other words: 


E(t) « t, (b —> ) (19) 


Thus, if we were able to determine £(t) for dispersion 
times from 0 to © we would expect & to vary with ¢ 
to the 2, 3 and 1 power, providing an inertial sub-range 
of eddy sizes exists and the rate of energy dissipation 
is constant. In the latter case the rate of change of £ 
is constant with time. 

These relationships for the relative mean dispersion 
of two particles were later extended by Batchelor (27) 
to include the case of a cloud of many particles. Gifford 
(23) further indicated that equations similar to (16) 
and (17) can also be shown to apply to the case of 
average cloud dispersion over a large number of identical 
trials and thus the equations should be applicable to 
relative diffusion of dye patches. 

The variances computed from the present experi- 
ments were assumed to be proportional to the true 
mean-square values. The fact that three-dimensional 
isotropy did not prevail in these experiments is evident 
from observed differences in the rates of vertical and 
horizontal diffusion. As discussed previously, it was not 
known for certain that the turbulence was truly aniso- 
tropic in the horizontal plane. The problem that arises 
in this respect is one of separating the effects of dis- 
persion due to wind action and vertical velocity varia- 
tion from that due to turbulent diffusion. If we can 


assume these effects to be manifested in longitudinal 
dispersion only, it is possible to check the validity of 
Batchelor’s predictions by analysis of the lateral vari- 
ance as a function of time. 

Figure 12 shows a typical relationship between the 
variance and diffusion time. It is seen that oy? is pro- 
portional to the time to a power less than 2. The trend 
of the plotted points do not indicate any apparent 
tendency toward an increase in slope at the larger values 
of time. It is possible that the experiment was not 
carried out long enough to show such a trend. However, 
an estimate of the time t, as given by equation (18) 
can aid in the analysis. Assuming the initial radius of 
the dye patch (= 6 ft.) to be proportional to the initial 
particle-separation in Batchelor’s predictions, and that 
a reasonable estimate of the constant rate of energy 
dissipation in the ocean is from 0.002 to 0.08 cm?/sec® 
(Okubo, (9) ) we have: 


iy (r)2/8 


a (oie = 4.3 to 1.3 minutes 


ti 


If these estimates are at all accurate, we can assume 
the time period during which (16) and (16a) apply is 
very small. Therefore, cy? would be expected to vary 
with diffusion time to a power ranging from 3 to 1. 

These data appear to indicate that the extent of the 
inertial subrange was limited to very small eddy scales 
in this particular oceanic area. It is possible that the 
state of turbulence in the area around the Orange 
County outfall is not ‘fully-developed” due to such 
complicating factors as the proximity of the shoreline, 
the shallow (50 ft.) depth, or the highly stable condition 
of the water column. Another important factor is the 
effect of the turbulent energy added by the presence of 
the boil. As Batchelor (26) pointed out, the range of 
eddy sizes over which the similarity hypotheses applies 
must be free of external influences, the basis behind the 
similarity theory being that a state of equilibrium of 
the smaller eddies must exist. This implies that the 
only source of energy for this range is that transferred 
from the larger eddies. If additional energy is fed into 
the range of small eddies, equilibrium may not obtain 
and the similarity theory may not apply. Even if the 
energy added by the effluent jet at the boil is of minor 
importance, it is still possible that the rate of energy 
dissipation within the field was not constant. As de- 
scribed earlier, the wind velocities varied considerably 
during individual test runs. Normally the wind speeds 
increased steadily from the early morning to late after- 
noon hours. The effect of a variable rate of energy 
dissipation on the rate of mean particle-separation can 
be seen by inspection of equations (16) and (17). 

The data summarized in Table 5 indicate that 
similar results were obtained for all tests conducted 
within the waste field. Additional ‘“‘control’’ experi- 
ments were conducted in coastal areas characterized 
by greater depth, and the absence of an effluent boil, 
and in most cases lower wind speeds. The results of 
these experiments are labeled c through f as shown in 
Table 5. The analysis of the data for cy? yielded results 
similar to those obtained for the experiments in the 
waste field. The empirical expressions for cy? show that 
the maximum exponent on ¢ was 2.36, all other values 
being less than 2. 
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10. DIFFUSION COEFFICIENTS 


It is possible to express the results of the diffusion 
analysis, at least in the x and y directions, in terms of 
apparent coefficients of diffusion. Thus, according to 


Okubo (6): 
d(sx.y? 
Ky = 3g Cee) (19) 


It can be shown from the equations given in the previ- 
ous section, that the so-called ‘4/3 law” of eddy diffu- 
sion is analogous to equation (17a), 1.e., differentiating 
(17) we have: 


d(const et’) 
i ee 


K= % ar = const et? (20) 


providing ¢ is not a function of ft. 
Eliminating t between (19) and (20) yields: 


K = const. (e)!/* (c)4/8 (21) 
Assuming the effective eddy scale | of the phenomenon 
to be proportional to o, we have 

K = const. (e)!/3 (1)4/8 (22) 


which, for a constant rate of energy dissipation e¢, is 
the ‘4/3 relationship.” 


When diffusion’ is controlled by larger and larger 
eddies than those found in the inertial subrange, the 
exponent on / decreases to zero, K becoming constant 
when the largest eddies characteristic of a given area 
control the diffusion process. The results of the dye 
patch studies appear to indicate that the ‘4/3 law” of 
eddy diffusion does not apply in the nearshore areas 
investigated in these experiments. 

Values of the average horizontal coefficients for each 
experiment were computed by use of equation (19), 
the empirical equations for the coordinate variances, 
and the average diffusion time. The corresponding 
average eddy scales were also computed by assuming 
1 = 4o,y. The computed values of K, ranged from 2.5 
to 1385 ft?/sec for average eddy scales of 650 and 3700 
ft. respectively. The corresponding values of Ky ranged 
from 0.2 to 7.4 ft?/sec for scales of 200 to 1000 ft. 

It is doubtful that the same procedure can be used 
to estimate the average vertical coefficients in the sur- 
face layers because of the influence of stability on verti- 
cal diffusion. The data in Table 5 shows that on several 
days o,2 was observed to be proportional to the diffusion 
time to a power less than one. The preceding discussion 
indicates that in this case K, would become smaller 
with time, which is contrary to theory. Regardless of 
these difficulties average values of K, were computed 
for the sake of comparison, and were found to range 
from 0.40 to 14.0 X 10-3 ft?/sec. 


B. DILUTION AND DISPERSAL UNDER CONDITIONS OF 
CONTINUOUS RELEASE: DYE PLUME STUDIES 


One of the main objectives of this research has been 
to study eddy diffusion of a surface waste field resulting 
from the continuous discharge of material from sub- 
marine outfalls. The equations for diffusion from a 
continuous fixed source can be obtained by integration 
(with respect to distance) of the equations for relative 
diffusion of a patch of material. Thus, information 
gathered on the horizontal and vertical diffusion of a 
dye patch in terms of the variances of the concentration 
distribution, can be used to predict or estimate the 
material distribution given by so-called steady plume 
models. However, this approach is limited (Okubo (9) ). 
Such complicating factors as the effects of local bound- 
aries and the non-homogeneous, non-stationary and 
anisotropic properties of a particular turbulent field 
must be given consideration. 

For these and other reasons it was considered ad- 
visable to conduct a number of field studies involving 
the continuous release of tracer dye from a fixed point. 
The studies were all conducted within the waste field 
at Orange County. 


1. FIELD METHODS AND SAMPLING 
TECHNIQUES 


The usual procedure followed in these experiments 
was to first pre-mix a known quantity of rhodamine-B 
dye with methanol for specific gravity adjustment. The 
dye solution was mixed in a 55-gallon drum mounted 
in a small skiff. Sea water was added to produce ap- 
proximately 45 gallons of final tracer solution. A small, 
calibrated positive-displacement pump with electric 
motor and variable speed drive was used to discharge 
dye at a constant rate. Power for the pump was fur- 
nished by a small portable generator. The skiff was 
anchored close to the boil by use of two bow and two 
stern anchors. Proper anchoring procedure was im- 
portant to minimize movement of the skiff and initial 
distortion of the dye plume. 

The typical sampling procedure employed during the 
continuous release experiments was similar to that used 
for the dye patch experiments with the following 
modifications: 


1. The dye pump was turned on and the discharge 

calibrated by volumetric measurement. 

It was assumed that steady state conditions were 

established within one to two hours after the ex- 

periment was begun. 

. The sampler was set at the depth of maximum dye 
concentration and repeated passes were made up 
the central axis of the dye plume. 

. Sampling across the plume axis was accomplished 
by making successive passes at known distances 
downstream from the skiff. Distances were esti- 
mated by use of a calibrated log-line velocity 
indicator and stop-watch. 


2. 


( 45 ) 


5. Vertical samples were taken periodically at known 
distances along the center line of the plume. 


Figures 26, 27, and 28 show a series of composite 
photographs taken during the experiments conducted 
on 27 and 28 June 1963. These Figures show the waste 
and dye fields as they appeared during two typical 
continuous release experiments and will be continually 
referred to in the discussion which follows. 


2. MATHEMATICAL MODELS 


Several statistical models have been proposed which 
describe the distribution of material in a steady state 
plume emanating from a fixed point source. The ma- 
jority of these models are based on the assumption that 
the turbulence field is isotropic, homogeneous, and 
stationary, and that the plume model is assembled by 
the superposition of an infinite number of overlapping 
Gaussian patches of material as shown in Figure 29a. 
If diffusion in the direction of the mean current is 
neglected, the steady plume can then be considered to 
be of the “disk-element’’ type proposed by Frenkiel 
(30), and shown in Figure 29b. 

Gifford (31) pointed out that a real plume presents 
a more complicated appearance than that shown in 
Figure 29. He indicates that the motion of the individual 
plume elements consists of two parts: spreading within 
individual elements, and meandering of the center of 
the element. Spreading is thought to result from turbu- 
lent eddies of sizes equal to or smaller than the plume 
element itself. The meandering is caused by large scale 
eddies which have no direct effect on diffusion taking 
place within the patch. Figure 30 depicts the ‘‘fluctuat- 
ing plume” models as visualized by Gifford. The prob- 
lem involved in attempting to predict the mean material 
concentration at any point with respect to the fixed 
axes, is obvious. Gifford does, however, propose a series 
of two-dimensional models that take into account the 
effects of both spreading and meandering of the plume. 
These models are based on an assumed Gaussian dis- 
tribution within individual patch elements as well as 
an assumed Gaussian form for the frequency function 
associated with the variability of the cloud center 
about the fixed axis. The instantaneous plume con- 
centration is given by: 


C(x,y,2) 


= gay oP [2 me | os 


where Q is the continuous steady rate of discharge of 
material from a point source, U is the mean current 
velocity, and ¢? is the average variance of the material 
distribution within the individual disk elements. e? as 
a function of dispersion time (or distance) was discussed 
in a previous section on relative diffusion. The disk 
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POINT SOURCE 


FIGURE 29a. STEADY PLUME FORMED BY SERIES 
OF OVERLAPPING "GAUSSIAN" PATCHES 


FIGURE 29b. STEADY PLUME FORMED BY SUPERPOSITION 
OF OVERLAPPING "GAUSSIAN" DISK ELEMENTS 
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POINT 
SOURCE 


FIGURE 30a "FLUCTUATING PLUME" MODEL COMPOSED 
OF OVERLAPPING "GAUSSIAN" PATCHES 


FIGURE 320b “FLUCTUATING PLUME" MODEL COMPOSED 
OF OVERLAPPING "GAUSSIAN DISK- ELEMENTS 
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element assumption, as indicated above, amounts to 
neglecting diffusion in the direction of motion. The 
concentration C at any point and time is considered to 
be a random function because of the variability of the 
position of the center of the individual plume elements, 
i.e., variability in the terms D, and D,. The equation 
for the mean concentration at any point is given by: 


ep — Loe py] CH 


where r = (y? + z?)!/2 and D? is the average value of 
the variance of the frequency function for the variability 
of the center of the element. This model also implies 
that diffusion within the individual elements is isotropic 
(equal in both the y and z directions). 


Gifford also proposes a two-dimensional model for 
the case where diffusion takes place anisotropically, 
i.e., at different rates in the two coordinate directions. 
The equation for the instantaneous concentration at a 
point is: 


C(x,y,2) 


C(x,y,2) ae Ge Det 


exp — 


4 Q [eee Be + G52) 
On (ay27,2) *2U Day? 
(25) 


The mean value of the concentration at a point is given 
by 


; Q 
CO) = te? + DAG? + DAU OP 


y? Ag 

ie + D,?) aa 2(6,2 + DD 20) 

The sampling procedures employed in these experi- 
ments involved measurement of the dye concentration 
distribution along the center line of the dye plume. 
The position of the center line was observed to change 
over long periods of time. However, the variation was 
small during individual sampling runs (see Figures 27 
and 28) and for this reason the results of these experi- 
ments have been interpreted on the basis of one- and 
two-dimensional steady plume models of the “disk 
element type’’. 


Equation (25) can also be used to describe the average 
concentration at a point in a steady state plume diffus- 
ing at the ocean surface. Using the method of reflection 
at the ocean surface, it can be shown that an approxi- 
mate equation for the mean concentration at any point 
in the plume is given by 


= 2 
C(xy,2) + er) 8 


Where the variances are a function of diffusion time or 
distance from the source, i.e., by making the substitu- 
tion t = x/U, the variances can be expressed approx- 
imately as functions of distance. 


Re Q y 
7 1 (oy2a,2 pee exp — (se 


The maximum theoretical concentration as a function 
of distance x can be obtained from (27) by setting 
Vy = 2 = zero and: 


Cmax = (28) 


ee 
1 (aye,7) 1/2Uy 


The one-dimensional form of (28) is obtained by 
assuming that diffusion takes place in the lateral 
direction only: 


d 
Cmax = a ay (29) 


where d is the average depth of the diffusing substance. 
Equation (29) is essentially the same as that proposed 
by Batchelor (32) for one-dimensional diffusion of 
material in the atmosphere. 

According to Okubo (9) it is possible that the average 
spatial distribution of diffusing substance might not 
be Gaussian, as assumed in the preceding equations. 
He indicates that Schénfeld’s (10) general equation of 
the superposable type may be applicable in this case. 
For one-dimensional horizontal diffusion from a con- 
tinuous fixed source, Schonfeld proposes the following 
solution: 


“a OU w(x/U) 
C(x,y) = i eaea (ett ees (30) 


where ‘‘w” is the mean diffusion velocity (assumed 
constant) and where the dispersion time ¢ has been 
replaced by its equivalent x/U. The maximum theo- 
retical concentration along the center line of the steady 
plume is obtained by setting y = zero in (30), 


eee (31) 


TWX 





The above discussion is not intended to be all- 
inclusive and many mathematical models proposed for 
diffusion under isotropic and anisotropic conditions 
(Cf. Sutton (83) and Fleishman and Frenkiel (34) ) 
have not been covered here. However, the models 
presented appear to provide the degree of flexibility 
and accuracy required to describe the concentration 
distribution of a quasi-steady state dye plume. 


3. EXPERIMENTAL RESULTS 


The data from the experiments conducted on 10 
November 1963 and 8 February 1964 at Orange County 
will be used to test the applicability of the above 
models. The data for these and other continuous release 
studies are summarized in Tables 6 and 7. On 10 No- 
vember approximately 6 pounds of dye per hour was 
discharged from the skiff anchored about 500 feet down- 
stream from the boil. The resulting dye plume resembled 
those shown in Figures 27 and 28. The dye and waste 
fields moved approximately parallel to the coast line 
in a southeasterly direction as shown in Figure 31. The 
average current velocity obtained by periodic radar fixes 
on the leading edge of the plume (corrected for diffusion 
effects) was 0.27 knots. The average wind speed re- 
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corded at the Southern California Edison Shore Station, 
was 1.6 knots. The average stability (S’ = 1230) for the 
0 to 3 meter depth indicated that very stable conditions 
prevailed in the water column throughout the test run. 
Vertical sampling of dye concentration along the center 
line of the plume revealed that vertical dispersion was 
small, the result of the high water stability and low 
wind speed. 

Figures 32 and 33 are log-log plots showing the com- 
puted values of the lateral and vertical variances of the 
dye concentration distribution as a function of distance 
from the source. They indicate that diffusion was es- 
sentially restricted to the horizontal plane. The maxi- 
mum dye concentration along the center line of the 
plume was determined for average times of 1.4, 2.2 and 
3.2 hours after beginning of the test. These values of 
Cmax were plotted as a function of distance from the 
source on log-log paper (Figure 34). The points indicate 
a linear relationship between the variables for distances 
of 500 to 5000 feet from the source. The least-squares 
trend line and corresponding empirical equation are 
also indicated. The empirical expression for Cymax iS 
assumed to be a reasonable approximation of the aver- 
age value of the maximum concentration as a function 
of distance. 

Equation (28), in terms of the variances of the con- 
centration distribution, was solved and compared to 
the observed values of Cmax. Thus, from the empirical 
expressions indicated in Figures 26 and 27, we have 
from (28): 


Cae 6.1 


165(x)—!-°! (ppm) 


(lb/ft*) 


The above equation is also shown for comparison in 
Figure 34. The estimated values of Cmax from the two- 
dimensional model agree reasonably well with the 
observed values. Since vertical diffusion of the tracer 
was minor, equations (29) and (31) were also used to 
estimate Cmax. Applying (29) with an estimated depth 
of the dye layer at 1.5 feet yields 


6.1 
2.51 [0.032(x)?-® ]!/2(1642) (1.5) 


53 .4(10)—4(x)—°-9°5 (1b /ft?) 
83 .5(x)—°-995(ppm) 


Cau = 


Using an estimated value for the mean ‘diffusion 
velocity” of 1.7 ft/min, as discussed in the previous 
section on relative diffusion, equation (31) gives 


6.1/60 
3.14(1.7)(1.5)x 


0.0127 (x)—!(lb/ft?) 
198(x)—!(ppm) 


The results of the above analysis are also plotted in 
Figure 34. It is apparent that for this particular set of 
data better agreement with the observed values of Cmax 
was obtained using the one-dimensional equations. 


ay ea 


I 


The importance of the vertical component of diffusion 
on the rate of change of Cmax can be seen by analysis of 
the data collected on 8 February 1964. The mean 
current velocity for this particular experiment was 0.14 
knots, considerably lower than that for the previous 
experiment. The mean wind velocity was 5.6 knots and 
the average value of stability, S’, was 375. For these 
conditions the extent of vertical mixing (Figure 35) 
was significantly greater than that found on 10 No- 
vember. 

The results of the analysis using the one- and two- 
dimensional models is shown in Figure 36. In this case 
the values of Cmax computed by use of (28) gave better 
results than those obtained from the one-dimensional 
models. Results similar to those discussed above were 
obtained in the other continuous release experiments. 


The results of these experiments point out that 
vertical diffusion becomes significant on days of rela- 
tively high wind speeds and low water column stability. 
This finding is in agreement with the results obtained 
in the studies of relative diffusion. 

It might be well at this time to discuss the assump- 
tions which form the basis for the above diffusion 
models and indicate how these conditions were satisfied 
in the present experiments. All of the models discussed 
above are valid for the case of diffusion from a contin- 
uous steady point source, in a field of homogeneous 
and stationary turbulence. Equations (28) and (29) are 
assembled by superposition of an infinite number of 
“Gaussian disk elements’’, in which diffusion in the 
direction of motion is assumed negligible and the ma- 
terial distribution within the element is assumed 
Gaussian. Equation (28) further implies that the diffu- 
sion is anisotropic, that is, diffusing at different rates 
in the y and z directions. Equation (31) includes a 
term for the ‘‘mean diffusion velocity” in the lateral 
direction which is assumed constant in space, thus pre- 
supposing homogeneity of the turbulent field. 

In these experiments the dye was discharged at a 
steady, continuous rate from a 1% inch diameter plastic 
tube located a few inches above the water surface. The 
initial width of the resulting dye plume was observed 
to be less than a foot. Thus, the requirements of steady 
continuous discharge from a point source were closely 
approximated in the field experiments. The require- 
ments of stationary and homogeneous turbulence, as 
discussed in a previous section of this report, were not 
necessarily satisfied in these studies. At times of high 
effluent discharge the water surface in the vicinity of 
the boil was extremely turbulent. Energy contributed 
by the rising effluent jet was added to that generated 
by local wind and wind wave action, and appeared to 
dissipate rapidly as the field moved away from the area 
of the boil. It is not known just how far this added 
energy is felt downstream, but field observations indi- 
cate it is relatively short-lived, on the order of a thou- 
sand feet. The rate at which energy was added by wind 
action was also variable. The wind velocity was usually 
low at the beginning of a test run, building up steadily 
during the day. Not much can be said about the trans- 
lational (homogeneous) properties of the turbulence 
fields encountered; further comment will have to wait 
until the turbulent properties in this, and similar 
oceanic areas, can be measured directly. 
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The assumption that diffusion within the dye plume 
is negligible in the direction of advective motion, is 
reportedly valid for conditions where the magnitude 
of the mean current is much greater than the turbulent 
velocity in the same direction. Ichiye (12) analyzed 
diffusion data collected under similar oceanographic 
conditions in Santa Monica Bay. He reports that a 
ratio of the mean current velocity to the turbulent 
velocity of approximately 10/1 was obtained. On the 
basis of Ichiye’s work it appeared reasonable to assume 
diffusion within the dye plume to be essentially one- or 
two-dimensional. 

The validity of the presupposition of a Gaussian dye 
concentration distribution was discussed in the section 
on relative diffusion. Plots of the instantaneous distri- 
bution of dye in the x, y and z directions for a typical 
patch were shown. Figures 37 and 38 are similar plots 
denoting, in this case, the mean values of the dye con- 
centration distribution observed while making repeated 
passes along the lateral and vertical axis of the plume. 
The figures also show the estimated value of the mean 
variances as computed from the average concentration 
distribution, as well as the hypothetical ‘‘Gaussian”’ 
distributions computed from the variances and Cyax. 
These data again point out that the Gaussian assump- 
tion is at least a good working rule for computations 
of this type. 

The agreement between the computed and observed 
results as shown in Figures 34 and 36 can be considered 
an indication of how well the required conditions were 
satisfied. 


4. MASS BALANCE ANALYSIS 


A mass balance analysis for all continuous release 
studies was accomplished by solving for Q (or Q’) in 
equation (28). Values of Cmax and the diffusion para- 
meters were computed from the empirical equations 
derived from the field data. A value of x corresponding 
to the average value over the test run was used in all 
such computations. The data summary (Table 7) shows 
that values of Q’ ranged from 144 to 294 per cent. 
These data agree with those obtained for the dye patch 
studies conducted within the waste field, and appear to 
indicate the influence of stability on the computed 
values of M’ and Q’. 


5. ADDITIONAL INVESTIGATIONS 


The results of the continuous release experiments 
indicate, to some extent, the applicability of equations 
(28), (29), and (81) to the case of diffusion from a 
quasi-point source. However, certain basic difficulties 
arise when attempting to apply these or similar models 
to diffusion of a large scale surface field in the sea. For 
example, point source models of the Gaussian type 
break down at small diffusion times (or distances), that 


is 
a —»> 0, and Cmax(0) —> ~ as t (or x) —>0 


The variances in the “disk-element”’ plume models 
can be interpreted in light of Taylor’s (35) theory of 
diffusion by continuous movements (one-particle type 
analysis). According to this theory: 

a « t,ast—}> o (32) 
g*cet7 fort <<a (33) 


where L; is the time scale of the turbulence (see Refer- 
ence (35) ). According to Gifford (31), the variance in 
Gaussian ‘‘disk element” steady plume models can also 
be interpreted in accordance with Batchelor’s predic- 
tions for relative diffusion as discussed earlier. Thus, 
we would expect the variance to be proportional to the 
dispersion time to the 3, 2, and 1 power as diffusion is 
affected by larger and larger eddy sizes. Also from 
Batchelor’s predictions we would expect that the vari- 
ance at ¢t (or x) = zero to be equal to the variance of 
the initial patch. For these conditions Cmax(0) must 
have a finite value as long as the source is not a true 
infinite point. The fact that Cmax(0) was not infinite 
in these experiments can be shown, for example, by 
estimating the value of the initial concentration from 
the data collected on 10 November. Assuming the dye 
plume to have an initial width and depth of 1.0 and 
0.25 feet respectively, we have from equations (28) 
and (16a) 





‘ue 
Cn) = S@*O)x20) )7U 


6.1 
~ 3.14(0.94 X 10-%)(1.84 X 103) 


1.12(Ib/ft?) 


lt 


The concentration of dye in the solution tank on the 
skiff was approximately 17 lb/ft?. This indicates that 
an immediate dilution of approximately 15/1 would 
have to take place to produce the concentration com- 
puted above. 


An inspection of the data shown in Figure 32 indicates 
that for this particular experiment, the lateral variance 
was proportional (empirically) to the distance from the 
source to the 1.81 power. There was no evidence of an 
(x)® regime in this or any other continuous release ex- 
periment. This can be seen by inspection of the equa- 
tions for a? shown in Table 7. It is possible, as 
mentioned in the section on relative diffusion, that the 
experiments were not carried out long enough to dem- 
onstrate an increase in slope of the gy? (x) relationship. 
In this respect results of the continuous release experi- 
ments verify those obtained from the dye patch studies, 
and indicate that Batchelor’s similarity predictions do 
not apply in this area. Analysis of the data from aerial 
photographs for the full-scale waste field (see next 
section) also appear to substantiate this conclusion. 


Again, it is not possible to analyze the data for vertical 
diffusion in the above manner because of the compli- 
cating effects of a vertical density gradient. 


Photographs shown in Figures 27 and 28 indicate 
that near the boil the shape of the waste field was 
essentially parabolic. This implies thata,? « x, as given 
earlier by equation (32) for the case of large diffusion 
times. Assuming the waste material to be Gaussian in 
its distribution (as will be discussed later), it was pos- 
sible to estimate the a,” (x) relationship from the photo- 
graphs for 27 June. The analysis revealed that z,? 
« (x)!°3, for distances to 6000 feet from the boil. A 
plot of a? — ay? (0) indicated little change in the ex- 
ponent on z for distances greater than about 100 feet. 
These results, although approximate, also indicated that 
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neither equation (16a) nor (17a) was applicable to the 
full-scale waste field. 

In order to investigate the details of the diffusion 
phenomenon in the waste field as it actually occurs, an 
experiment was performed in which the entire waste 
effluent was tagged with rhodamine B. The dye was 
injected into the effluent at the Orange County Sanita- 
tion District’s Treatment Plant No. 2 prior to discharge 
through the outfall line. Since the flow from the plant 
is variable with time, it was necessary to proportion the 
rate of dye injection by use of a variable speed positive 
displacement pump. A constant dye concentration of 
2 ppm was maintained in this manner for approximately 
5 hours. During this period the flow rate from the plant 
varied from approximately 42 to 100 million gallons 
per day, the greatest increase in flow rate occurring 
during the last two hours of the experiment. 

Figure 39 is a composite photograph (oblique) taken 
approximately 1.5 hours after the dye appeared at the 
boil. The wind moved across the dye plume throughout 
the test and apparently contributed to the “layering 
effect”? shown. A visual estimate of the scale of the 
figure can be obtained by noting the position of the 
RV VELERO IV (106 ft. long) in the upper left-hand 
corner of the photograph. 

The complicated nature of the dispersion phenomenon 
for this particular day can be seen in Figure 40. This 
figure is a composite of several vertical photographs 
taken 3.5 hours after the start of the test. The overall 
complexity of the dispersion pattern for this experiment 
is further illustrated by comparison of the waste plume 
for this study with those shown in Figures 27 through 
28. The waste field was moving in approximately the 
same direction in all three experiments (see Figure 31). 
The reasons for the different flow patterns are difficult 
to determine. The magnitude and direction of the wind 
and surface currents were approximately the same (see 
Table 6). The rates of waste discharge, while variable 
during each day, were similar, since all three tests were 
conducted at approximately the same time. It was 
obvious that any attempt to describe the space-time 
distribution of dye for this particular experiment would 
be extremely difficult. 

The sampling procedure followed for this study was 
similar to that employed for the other continuous re- 
lease experiments. Samples along the longitudinal axis 
of the plume were made at various average time periods. 
Figure 41 shows the results obtained when Cmax was 
plotted against distance on log-log paper. The results 
show that the logarithm of the concentration drops 
gradually with distance from the boil and becomes ap- 
proximately linear for distances greater than about 1000 
feet. The least squares equation for the linear portion 
of the curve indicates that Cmax « (x)~°-*5, quite differ- 
ent than might be expected from the results of the 
previous experiments, as will be discussed later. 

The vertical distribution of dye concentration was 
determined at known distances along the visual center 
line of the plume, and the results are indicated in 
Figure 42. Plots of dye concentration with depth 
indicate the distribution was approximately Gaussian. 

One especially interesting aspect of this experiment 
was the lateral distribution of dye concentration. Pre- 
vious (visual) observations of the waste plume near 
the boil indicated that a sharp concentration gradient 


existed between waste and background water in this 
area. The effect appeared to extend downstream for 
several thousand feet. This observation was confirmed 
by dye concentration measurements as shown in Figure 
43. The plot indicates that there was a sharp increase 
in dye concentration as sampling proceeded across and 
toward the center of the plume at the boil. The concen- 
tration then leveled off and subsequently increased to 
a peak value at the center of the boil. A similar pattern 
was manifest on the other side of the plume. Figure 43 
shows that a distance of 1350 feet downstream the 
concentration distribution was quite different from that 
found at the boil. The fact that lateral dye concentra- 
tion distribution was non-Gaussian is also evident. 


The sharp drop-off in dye concentration at the lateral 
limits of the plume persisted for several thousand feet 
downstream and was somewhat more pronounced on 
the seaward side. The reasons for this effect are not 
known, but the phenomenon also has been observed by 
Gunnerson (18) around the old City of Los Angeles 
outfall in Santa Monica Bay. Gunnerson drew an 
analogy between this effect and those of ‘‘mixing’’ and 
“stirring,” originally proposed by Eckart (36). Eckart 
distinguished between “stirring” and ‘mixing’ pro- 
cesses, by considering “‘stirring’”’ as the distortion of 
fluid elements by a shearing current. Large concentra- 
tion gradients are produced by this mechanism, which 
promote “mixing” in the direction of the gradient. 
Supposedly these two processes interact to produce the 
overall observed diffusion effects. Further discussion on 
this phenomenon will not be given here except to state 
that the diffuser system at the terminus of the Orange 
County outfall line is relatively inefficient and probably 
is responsible for the sharp concentration gradients 
shown in Figure 43. The diffuser consists of a series of 
small circular openings on both sides of the effluent 
line near the outlet. The end of the line is open, and 
observations by divers indicate that approximately 
70% of the waste discharges through the open end. The 
result is poor initial dilution as the effluent jet rises to 
the surface. An outfall located in deeper water with a 
more efficient diffuser would provide better initial 
dilution, reducing the sharp concentration gradient at 
the effluent—sea water interface. It is possible that 
under these conditions the lateral concentration distri- 
bution within the waste field would be more nearly 
Gaussian. 


6. DIFFUSION FROM A VOLUME SOURCE 


The plot of Cuax (Figure 41) suggests that equations 
of the point source type could not accurately describe 
the process as shown for all distances from the source. 
Gifford (37) proposed a series of statistical models for 
the case of atmospheric diffusion from volume sources, 
which indicate a relationship similar to that shown in 
Figure 41. His proposed models are for instantaneous 
volume sources in an isotropic, homogeneous, and sta- 
tionary atmosphere and thus are not directly applicable 
to the case at hand. We will present a brief review of 
his work and try to illustrate how at least one of Gif- 
ford’s models can be extended to describe the results 
of the experiment of 17 June. 

Gifford assumed a volume source to be the sum of 
an infinite number of instantaneous point sources each 
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CONTINUOUS DISCHARGE 


ORANGE COUNTY OUTFALL 


17 JUNE 1964 


3.5 HOURS 


0.32 KNOTS 


COAST LINE 
APPROX. |.2 MILE 
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of strength Q (x’, y’, z’), located at points (x’, y’, 2’) 
in space. The mean concentration distribution of the 
volume source is found by integration with respect to 
x’, y’; and z’. Thus: 


C(x,y,z,t) = (2a Si Bie nh 4 Q(x oY. 2’) 
as \2 as \2 yee ‘\2 
Reserpine Nat by sb ob 8 2’)? | 
dx’dy’dz’ 
where @? is the variance which is a function of diffusion 


time (or distance) and which is assumed to be equal 
in all directions. Gifford further shows that for a 
volume source with an initial Gaussian concentration 
distribution, the above integration yields: 


C(x, y,z,t) (34) 


a M | 
~ Tr(Qe? + b=) 22 OXP are +b 


where C is the average concentration, r2 = x? + y? + 
z?, and 1/b is the variance of the initial Gaussian cloud. 
He points out that the result is identical, except for 
1/b, with other Gaussian point source models. Equation 
(34) can be looked upon as the case where an initial 
Gaussian cloud, with variance 1/b, is formed by some 
fictitious point source acting upstream (originally 
reported by Holland (38) ). 


Gifford did not extend his model to the case of the 
continuous release of material from a volume source. 
However, it seems possible to do so by analogy with 
equations (28) and (29). Thus, for a continuous release 
of material from a volume source, neglecting diffusion 
in the x direction (disk-element assumption) we have: 


C(x,y,z) 


na 2Q 
aU [2ay2 + ay2(0) ]!/2 [2¢,2 + 2,20) }!/? 


- Las taw + ETOH 2 


where the variances are now a function of x by sub- 
stitution of x/U for t. The maximum concentration 
along the center line of the plume is obtained by setting 
y = z = zero in (35). Thus: 


Oly — 


7Us,(0)z.(0)| Foy ae + ay iE sy +1]'"@6) 


The initial maximum central concentration where 
Xx = zero is: 


ie koleaecre bee (37) 


1 Uay(0)az(0) 


Equation (36) can be rewritten in terms of Co as: 


Cmax = aa i a i (38) 
aon ay?(0) ul [ 7,7(0) Ut 
The one-dimensional forms of equation (37) and (38), 


for diffusion in the horizontal plane, can be expressed 
by: 





Gini = =a (39) 
Exo (0) + 1 
and: 
De doshas | 
© = Ty9,(0) Se 


Where Cp» is again the initial central concentration in 
weight (or numbers) per unit volume, Q is the steady 
discharge rate in weight (or numbers) per unit time, 
d is the average depth of the waste field, and a, (0) is 
the standard deviation of the initial lateral waste 
concentration distribution. 


Application of the above equations to a typical design 
situation involves prior knowledge, estimated or mea- 
sured in the field, of the variance at the point of waste 
discharge, the lateral and vertical variances of the field 
as a function of distance from the origin, and the initial 
waste (or waste constituent) concentration. 


The data collected on 17 June can be used to (roughly) 
indicate the validity of the above equations. The 
problems involved in such an analysis are obvious, the 
discharge was variable and the lateral dye distribution 
was more nearly uniform than Gaussian. The data 
indicate that the lateral spread of dye was minimal; if 
anything, the plume decreased in width with i increasing 
distance from the boil. However, the vertical spread of 
the waste increased with distance as indicated in 
Figure 42. These conditions suggest that the one- 
dimensional equations, with the average width W re- 
placing d and a,” replacing ay”, could be applied. The 
estimated physical parameters for the test are: 


42 lb/hr = average for test run 
0.32 knot = 1950 ft/hr 
1500 ft. = average plume width with 
uniform concentration 
distribution 
1 ft. (estimated initial depth of the field 
is 2 ft and distribution assumed to 
be Gaussian.) 


4 GSI 
tou i 


a, (0) 


The initial central concentration, assuming reflection 
of the boundary for the vertical case, is: 


GO = y= 
°  Uvxa,(0) 


16 .3(10)®(1b/ft?) 
254(ppb) 


This value should be compared to a measured value of 
approximately 265 ppb at the boil. The maximum 
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concentration along the center line of the waste plume 
is: 
Co 

26,7 1/2 
[ 300) ial 
where o,” is given empirically in Figure 42. The results 
of this analysis are shown in Figure 41. The fit of the 
theoretical curve to the observed concentration distri- 
bution is good. However, the results of this analysis 
obviously can not be taken as confirmation of the 
proposed equations, and additional data on large scale 
diffusion of surface waste fields is needed for further 
analysis. 

In most cases the horizontal component of diffusion 
within the waste plume would be the most significant, 
if not the controlling, factor in the overall dilution 
process. Therefore, the volume source equations used 
to predict Cmax would be, including vertical diffusion, 
of the two dimensional type. An example of the method 
of analysis applicable to this case will aid in the present 
discussion. Using the data collected for the experiments 
of 27 June 1963 and 17 June 1964 as a guide, we can 
solve the following hypothetical problem. 


Com = 


Basic assumptions: 

1. Diffusion of the waste plume is two-dimensional, 
ie., in the Y and Z directions. Longitudinal 
diffusion is assumed to be negligible. 

. The diffusing substance (dye) is conservative. 

. Diffusion will be at or near the ocean surface. 

. Steady state conditions exist. 

. The lateral and vertical distribution of diffusing 
material is Gaussian in form. 


Ot GC bo 


Physical Conditions: 
2 = 42 lb/hr (dye) 
= 2200 ft/hr = 0.36 knots 
=~ 45 ft = width/4 
7,(0) += 1ft = depth/2 
o771=1160(x)j ft? 
Fe 7 GH 10 Fi x) 9 rit 


Applying equations (37) and (38) gives the following 
results: 


Co = 2.8(ppm) 


2ey* = 26,” 1/2 
Exo y?(0) +1] ae 2° (0) +1] 
Substitution of the assumed empirical equations for the 
coordinate variances into the above equations yields 
the relationship for Cmax. The results are shown in 
Figure 44 and are labelled “example A.” 

The point source equation for Cmax for the assumed 
conditions can also be solved for comparison. Thus, 
substitution of assumed values for the variances into 
equation (28) gives 


Cmax = (ppm) 


42 
~ (2200)r [60 X 7.6 X 10-3]1/2(x)0-%7 


0.903(10)—2(x)—°-9"(Ib/ft?) 
141.1(x)-°-*"(ppm) 


Aree (x) 


I 


The results are also shown in Figure 44 as example B. 
A comparison of the volume and point source equations 
indicate that the equations converge, giving essentially 
the same results, for X > 1000 ft. However, for dis- 
tances between 0 and 1000 ft. the point source equation 
gives values of Cmax that are much too large. These 
results indicate that for diffusion near the source, the 
volume source equations give a more accurate and 
realistic estimate of waste dilution. 

The effect of neglecting the vertical component of 
diffusion in the above example can be seen by solution 
of equation (39) with Co = 2.8 ppm for purposes of 
comparison. Thus 


Ls) 1G 1] 


The resulting expression for Cmax(x) is shown for com- 
parison in Figure 44 as “example B.” The effect of 
neglecting vertical dispersion in this hypothetical ex- 
ample is obvious. 

The assumed expression for the lateral variance 
indicates that in this example the lateral spread of dye 
(or waste) is parabolic, as indicated in Figures 26, 27 
and 28. The spatial distribution of dye concentration 
can be determined by solving equation (35) for values 
of y (at given values of x) corresponding to any arbi- 
trary dye concentration. Figure 45 shows a series of 
iso-concentration contours computed by use of equation 
(35) and the values of Co and oy? given in example A 
above. 

Whether g,? is proportional to x to the first, second 
or third power will depend on the initial size of the 
volume source as compared with the range of eddy 
sizes (eddy spectrum) characteristic of a given oceanic 
area, aS suggested by Orlob (39). 

The proposed volume source equations are obviously 
approximate and need further verification. Additional 
full-scale data on horizontal and vertical diffusion from 
existing outfalls is needed. 


The proposed equations if valid, give the engineer 
an additional means of estimating the dilution and 
physical spread of a surface waste field. It should be 
pointed out that these equations offer a means of esti- 
mating the initial central dilution (1/Co) of the waste 
at the boil, as well as subsequent minimum dilution 
(1/Cmax) due to diffusion of a steady waste plume. The 
equations can be modified to include the diffusion of 
microorganisms (coliform bacteria for example) by 
inclusion of a term for bacterial disappearance (die-off, 
settling, etc.). Since the equations developed for Cmax 
are based on probability density relationships we can 
assume that for the case of microorganisms 


Cae = 





ear (41) 
— (XxX 
~ F Qay? Qa yaaa [Ka)] 
lam + 1] (aes 1 
Where the last term in (41) is the so-called ‘“die-off” 
factor as a function of diffusion time (t = x/U), and 


67 


FATE OF WASTES IN THE MARINE ENVIRONMENT 


000‘'0o! 


o000'o! O00S 


) adwoxg 


bp SYNSIS 


(13334) SJONVLSIC 


ooo! 


00S ool OS Ol G | 
10°O 


SISAIVNY S30YNOS SAWN IOA 


JOYNOS WOYS SJONVLSIC 
“SA 
NOILVYLNSONOD WNWIXVW 


y 9ejdwoxg 


(894N0S juI0d) 
gq 9j\dwoxg 





(NWdd) NOILVYLNADNOD 


FATE OF WASTES IN THE MARINE ENVIRONMENT 


68 


(,01 * 1334) 
6 


Sb 3yNDdIS 


JONVLSIGC 
S 


Wdd 8'2=9 


Gdd OS! 


< 
| 
7 
fn 
+ 
x 
Xe) 
ce 


@ddOol 


NOISNSSIIG IWNOISNSWIG—-OML 
JOYNOS SAWNIOA SNONNILNOD 


AWN Id SAC 





FATE OF WASTES IN THE MARINE ENVIRONMENT 69 


the concentration is expressed as number of micro- 
organisms per unit volume. 

The material presented in this section of this report 
is to aid in the development of accurate mathematical 
formulations that can be used in the efficient and 


economical design of future submarine outfalls. In this 
respect the material presented herein is intended to 
supplement the excellent previous work of Brooks (40), 
Pearson (42) and Rawn and Palmer (41) on this subject. 


C. PRIMARY PHYTOPLANKTON PRODUCTIVITY AND 
PRODUCTION IN A MARINE WASTE FIELD 


1. INTRODUCTION 


The utilization of natural bodies of water as recep- 
tacles for unwanted wastes predates the written history 
of man. With the growth of communities and industries 
the increase in the quantity of such wastes resulted in 
drastic changes in the quality of the receiving waters 
and concurrent inevitable changes in the biota. 

The biotic changes, first noted as the alteration, 
depletion, or disappearance of certain faunal popula- 
tions, are still most frequently assessed in terms of 
standing crops, i.e., the presence and quantity of various 
biotic elements. Examples of such observations include 
the increase in quantity of tube worms and the decrease 
of benthic algae as outlined by Bartsch (43), the ap- 
pearance, subsequent loss, and later recovery of a 
commercial fishery (Hubbs (45) ), and alterations in 
biological associations (Hartman (44), Tibby and 
Barnard (46) ). 

Understandably, a great proportion of such studies 
have dealt with the effects of sewage disposal in bodies 
of fresh water such as rivers, lakes and streams. Much 
of this work has been summarized by Hynes (47). 
Confining water masses and restricting currents to well- 
defined channels permits the build-up of concentrations 
of sewage effluent and its products to a quasi-steady 
state condition along the water course downstream from 
the outfall. As a function of distance and hence of 
dilution and alteration of the sewage, a series of inter- 
graded ecological niches are found in which populations 
characteristic of the changed conditions will supplant 
earlier communities. Observations of this nature have 
led to variously defined systems of biological ‘indicator 
species” to evaluate levels of pollution. Those systems 
are reviewed by Hynes (loc. cit.) and will not be dis- 
cussed in detail here other than to indicate that any 
single system may not be universally applicable because 
each stream, river or lake may have different normal 
populations and different ecological conditions. 

The differences between fresh water and marine 
environments require some modification of the methods 
used in their biological study. Some of the methods 
for marine work are detailed by Jannasch and Jones 
(48), Lackey (49) and Barnard and Jones (50). 

Studies in the southern California marine environ- 
ment of standing crops of organisms (Clendenning and 
North (51), Reish (52), Resig (53), Bandy, Ingle and 
Resig (54, 55), Hartman (44), Tibby and Barnard (46), 
Allan Hancock Foundation (56) ) have described ar- 
rangements of biotic communities in relation to outfalls, 
arrangements which are similar in principle to those 
in fresh water environments. 

Assessment of the influence of sewage in fresh and 
sea water on biological processes such as photosynthesis 
and respiration often has been made on the basis of 
standing crop or BOD. However, direct measurement 
of the processes involved, and particularly in algal pop- 


(70 ) 


ulations, has been utilized only in isolated instances 
and for specific purposes. Among such studies are those 
of Goldman (57) who considered the role of sewage in 
accelerating the natural eutrophication of Lake Tahoe, 
and Lackey and Sawyer (58) and Lackey (59) who 
concluded that nutrients added by sewage contribute 
to the growth and blooms of planktonic algae. Steeman- 
Nielsen (60) believed that the productivity of a Danish 
lake fertilized by nutrient salts of sewage origin was 
limited by damage to the plants caused by a rise of pH 
engendered by CO, depletion. However, Eppley and 
MaciasR (61, 62) have found that under certain condi- 
tions an organic source of carbon may be utilized in 
photosynthesis by sewage pond algae. Algal blooms in 
small embayments have been attributed to the nutrients 
added by polluted water (Ryther (63) ) although for 
open coastal waters there are apparently no factually 
documented reports. 

Standing crop measurements give a summation of 
the momentary states of balance between the processes 
that are occuring but do not give an accurate picture 
of the physiological effects of the sewage on the biota 
or any indication of the rate of turnover of the biota or 
of the organic and inorganic nutrients contributed by 
the sewage. Such information on rates are of first im- 
portance in considering the effective capacity of an area 
to handle sewage without profound alterations in the 
biota, and in predicting the effect of an outfall on the 
biological resources of the area. It has been the intent 
of this project to utilize existing methods and to develop 
new techniques in an attempt to assess this magnitude 
of biological and chemical turnover within a marine 
sewage field. 


2. AREAS STUDIED, AND METHODS 


The area selected for study lies off Newport Beach, 
California, where the Orange County Treatment Plant 
#2 discharges about 75 mgd of mixed domestic and in- 
dustrial wastes through a pipeline terminating about 
114 miles off the mouth of the Santa Ana River in 
about 60 feet of water (Figure 31). The sewage, most 
of which has had only primary treatment, is warmer 
and less saline than sea water and rises to form a visible 
boil on the surface. At the surface the water contains 
about 5-10% sewage effluent. Currents in the area 
consist largely of longshore currents running southward 
roughly parallel to the coast, although occasionally a 
northward flowing current is present. Frequently there 
is an onshore movement of the surface waters in addi- 
tion to the longshore currents. The general southward 
transport of surface water along the coast is disrupted 
about 3 miles south of the outfall line by turbulence 
associated with the Newport Submarine Canyon. 

This area was selected over others on the southern 
California coast for several reasons. First, the surface 
boil and the downstream sewage field were nearly always 
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easily seen, in contrast to the situation in Santa Monica 
Bay where the Hyperion effluent field is much larger 
but is extremely difficult to locate. Second, the type 
of mixed domestic and industrial effluent discharged is 
fairly representative of the other major outfalls along 
the coast. Third, the Orange County outfall was within 
convenient sailing time from the ship’s berth in San 
Pedro. 

Except for some special purposes, stations were occu- 
pied in relation to the trajectory of the diffusing sewage 
field. The method of choice in identifying and tracking 
specific water masses was to use the fluorescent dye, 
rhodamine-B, supplemented by drogues. Occasionally 
a dog-leg or grid pattern was established in relation to 
the boil and to the currents prevailing at the time. 

The use of current markers of one sort or another has 
the great advantage of permitting a study of the time- 
sequence of events along the trajectory of the diluting 
sewage field. Stations were occupied at timed intervals 
near the marker which (supposedly) stayed in the same 
patch of diluting sewage. The selection of appropriate 
types of current markers can be quite critical. Drogues, 
although easy to use and inexpensive, can give an in- 
accurate picture of the currents being measured since 
they are subject to wind drift and may be caught in 
eddies that are not following the main current. Dye 
patches or dye streams released at the boil are less 
subject to wind drift and, as they disperse to form a 
larger sized patch, are less subject to being retained in 
small eddies or other peculiar, small-scale features of 
the circulation. However, the patches ultimately diffuse 
beyond the level of easy detection, and they cannot be 
readily followed at night. The dye itself (rhodamine-B) 
has a slight adverse effect on carbon-14 uptake of 
phytoplankton, and for this reason productivity samples 
usually were taken just ahead of the dye field. The most 
effective current marker we have tried is to utilize the 
entire sewage field by continuous release of dye in the 
treatment plant prior to pumping it and the sewage 
effluent through the outfall line. In this case, the entire 
sewage field was highly visible and the relative concen- 
tration of sewage in the field was well defined. Because 
of the expense of a sufficient quantity of dye to provide 
a continuous-release tag adequate for extended periods 
of study, this method does not seem promising as a 
routine measure. 

For comparison with the study area, a control area 
along a 20 mile stretch of coast was selected off Camp 
Pendleton equidistant between Oceanside and San 
Clemente. Oceanside lies some 35 miles south of the 
study area. Although the submarine topography and 
general hydrographic conditions are not exactly like 
those of the study area, there was no other site within 
easy sailing distance that could be assumed to be free 
of sewage. Camp Pendleton, the nearest sizeable center 
of population, reprocesses its sewage and does not dis- 
charge any of it into the coastal waters. San Clemente, 
9 miles north and upstream is the nearest community 
utilizing the sea for waste disposal. The average volume 
flow of this outfall is 0.5 mgd., and its effect is appar- 
ently negligible in the control area. ; 

Stations in this control area were occupied either 
using a predetermined dog-leg pattern or lines of sta- 
tions at right angles to the shore. Sampling was carried 
out as time and boat scheduling permitted and yielded 


data representative of the situation to be expected in 
uncontaminated water along a stretch of coast roughly 
similar to the study area in size and hydrographic 
conditions. 

All samples for determining productivity and pig- 
ments were collected in non-metallic samplers to avoid 
the problem of poisoning by metals. For surface samples 
three gallon plastic buckets were used. Subsurface 
samples were collected using ‘snatch’ bottles, Van 
Dorn samplers, and a specially designed and constructed 
valve sampler. The “‘snatch” bottle was a weighted 
2-liter plastic bottle which was lowered on a hand line 
to the appropriate depth. A separate line connected to 
a rubber stopper in the mouth of the bottle was used 
to “snatch” the stopper from the bottle. The Van Dorn 
samplers with a capacity of about 6 liters were pur- 
chased from Kahl Scientific Instrument Co., San Diego, 
California. These samplers consisted of a plastic tube 
closed by rubber force cups held under tension by a 
piece of rubber tubing. The specially designed samplers 
were constructed to our specifications by Logan Smith 
and Associates, Harbor City, California, and were 
extensively used in the later stages of this study. They 
are described in the later section on chemistry. 

The method of assaying productivity used during 
this project was slightly modified from the isotopic 
carbon tracer method first described by Steeman- 
Nielsen (64). Each sample collected at each station was 
subdivided and a portion was used to completely fill 
duplicate series of clear and opaque 125 ml. glass- 
stoppered Pyrex bottles. These bottles were stored at 
ambient temperatures in a covered box. Addition of the 
isotopic tracer in a 34°/o) NaCl solution containing 2 
micro-curies of carbon-14 as sodium carbonate, and 
incubation was begun the following morning. Using a 
standard time of day for incubation eliminated the 
correction for diurnal periodicity in the photosynthesis 
of phytoplankton reported by Doty and Oguri (65). 
Incubation was carried out for production measure- 
ments by suspending the bottles at the depth and 
location of collection from a buoyed and anchored line, 
or for productivity measurements by placing the bottles 
in a box on deck where they were immersed in running 
sea water and exposed to ambient daylight. Following an 
incubation period of about 3-4 hours for samples incu- 
bated in the deck box, and 5-7 hours for samples 
incubated in situ, samples were filtered through 
Millipore AA filters and the filters were stored in a 
desiccator for return to the laboratory where the 
activity retained on the filtered samples was counted 
in a low-background flow counter. Counts from the 
clear bottles minus counts from the opaque bottle were 
used to compute milligrams of carbon fixed per day per 
cubic meter. 

Pigment analyses were carried out according to the 
method of Richards and Thompson (66). For this 
purpose, one-half to one liter of sample water was 
filtered through a Millipore AA filter. As the last of the 
water went through the filter, a small quantity of a 
suspension of MgCO; was added to retard breakdown 
of chlorophyll. The’filters were then stored in a darkened 
desiccator held at refrigerator temperatures. In the 
laboratory, pigments were extracted by placing the 
dried filters in 90% acetone for 12 to 24 hours in the 
refrigerator. The optical absorption at the wavelengths 
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of the peak absorption for the chlorophylls was then 
measured either in a Beckman model DU or DB spec- 
trophotometer with, respectively, 10 cm. or 4 cm. light 
paths through the sample. Pigment concentrations were 
calculated as mg. or MSPU of the different pigments 
per cubic meter according to the formulae listed by 
Richards and Thompson (loc. cit.). 


3. CHANGES IN PRODUCTIVITY 
ASSOCIATED WITH WASTE 


In the Annual Report for 1962-63 (1) the reasons for 
considering the productivity as a function of standing 
crop were discussed. Briefly, due to a persistent and 
patchy red tide in the study area the pattern of pro- 
ductivity found both in the background and the sewage 
field was variable and inconsistent. It was believed that 
a clearer and more readily interpreted pattern might 
be obtained if productivity were considered as a func- 
tion of the biological standing crop rather than as a 
function of volume or surface area of water. The 
inherent tedium of cell counts and species lists precluded 
their use in this survey. Photosynthetic pigments, the 
chlorophylls, were selected as a measure of the standing 
crop because of the ease and rapidity of their deter- 
mination and their more direct relationship to the 
productivity of the organisms. Although correlations 
between productivity and changes in chlorophyll con- 
tent in marine phytoplankton are not expected to be 
high, our work was confined to a relatively small area 
containing a population that can be assumed to be 
relatively uniform in make-up and in photosynthetic 
efficiency during any one cruise. 


Table 8 presents data from one cruise during which 
samples were collected both in the control area south 
of San Clemente, and in the study area off the Orange 
County Treatment Plant. Both areas are about the 
same relative size, extending 6—7 miles along the coast 
and 2-3 miles offshore. The data for all stations show 
correlation between productivity and total chlorophyll 


TABLE 8 


The correlation and regression of productivity per unit of chloro- 
phyll of surface and sub-surface samples at a series of stations in the 
control area south of San Clemente, and in the study area at the 
Orange County outfall, 31 October-2 November 1962. Stations in 
the study area are numbered consecutively in order of time occupied 
during drift from the boil with Station 1 located at the boil and 
Station 8 farthest from it. The time between stations was approxi- 
mately two hours. 








Study area 
Control 
area 
all All 


stations | stations Pairs of stations 











1&2 3&4 5 &6 7&8 


Station (0-2 (4-6 (8-10 | (12-14 

number 1-8 1-8 hr) hr) hr) hr) 
No. of samples____-_ 24 23 6 5 6 6 
Correlation_______- 0.925¢ | 0.581¢ |—0.625 0.089 | 0.890* | 0.9187 
Regression _______- 25.80 45.73 |—8.27 3.03 49.91 91.22 
* Significant 


+ Highly significant 


of 0.925, n = 24, at the control site and 0.581, n = 23, 
in the diluting sewage field at the Orange County area. 
The data for the Orange County area are further broken 
down to cover pairs of consecutive stations starting 
at the boil and ending 14 hours after discharge. The 
correlations then run from negative nearest the boil 
to high and positive away from the boil. The correlation 
at the pair of stations 7 and 8 farthest from the boil is 
similar to that in the San Clemente area. The lack of 
good correlation close to the boil reflects the confused 
picture and rapid change of the relationship in the first 
few hours following discharge of the sewage. 

The regressions of productivity per unit of chloro- 
phyll indicate that even at the last pair of stations in 
the field there was an enhancement of productivity 
compared to the control area while the correlations 
suggest that the sewage remaining no longer exerted 
the inhibiting effect that was prominent during the first 
six hours (Stations 1 & 2, 3 & 4). 

Figure 46, reproduced from the Annual Report for 
1962-63 (1), shows the production and productivity 
per unit of chlorophyll found in the surface sewage field 
during eleven cruises in the study area at a time of 
intense red tide. In this and subsequent Figures, and 
in the discussion, one may read ‘‘production” for 
“productivity”, where “production” is measured by 
carbon-14 uptake in seeded bottles suspended at the 
depths from which the samples were originally taken, 
and ‘productivity’ is the potential production deter- 
mined from samples incubated in the box on deck at 
ambient conditions of light and temperature. During 
the time-interval covered by each of the field studies, 
the sewage field remained at the surface and produc- 
tivity therefore closely approximated production in situ 
at the sea surface. A number of early experiments con- 
firmed the fact that productivity of surface samples 
and production in situ at the surface yielded almost 
identical results. This was a fortunate circumstance 
because time-consuming in situ measurements then 
became unnecessary for purposes of this study. 

Figure 46 shows a fairly clear pattern of low produc- 
tivity near the boil, an increase to a high some distance 
away from the boil, and then a dropoff towards back- 
ground values. The situation near Newport Canyon is 
somewhat confused because of interruption of the 
normal longshore currents toward the southeast. Al- 
though the intensity of the red tide in the area during 
this period varied, as long as it did persist the phyto- 
plankton of the area were essentially a unialgal popula- 
tion. The overall photosynthetic efficiency of the 
phytoplankton, therefore, did not appear to vary in 
time as it might under more normal conditions when 
populations of different relative composition would 
succeed one another. Consequently, subsequent cruise 
data showed a smaller and more variable rate of carbon 
fixation per unit of chlorophyll, although the same 
general sequence of low-high-background productivity 
occurred with remarkable regularity. 

The data for the later cruises (after May, 1963) when 
there was little or no red tide were grouped according 
to season and the resulting patterns showed some con- 
sistency. These seasonal patterns are presented in 
Figures 47, 48, 49 and 50. 

Neither the seasonal coverage nor the station patterns 
were as extensive in the control area as those in the 
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study area, but the data for the control area can be 
extrapolated to cover the same size and shape of area 
for purposes of comparison. The data on surface pro- 
ductivity (production) shown in Figures 47 through 50 
are tabulated by seasons in Table 9 together with 
similar data from the control area. In the study area 
it is evident that there was a seasonal pattern showing 
high production in the spring and summer and low 
production in fall and winter. It is of interest to note 
that a red tide in the study area added a second and 
greater maximum of productivity in the fall with pro- 
duction about 5 times higher than in non-red tide 
periods. In the winter and spring the effect of the red 
tide on total production was considerably smaller. Red 
tide productivity in the winter months was 65% higher 
than in the absence of red tide and in spring it was 23% 
higher. With or without red tide, winter production 
was relatively low. No red tide was evident during the 
summer cruises. 

The control area was investigated only in the spring 
and fall. For these two seasons the results were similar 
to those from the study area in that spring production 
was about double that in the fall. Red tide in the control 
area was evident in only a few isolated samples during 
the few cruises to that area and these samples were 
omitted from the calculated data in Table 9. 

A comparison of the non-red tide data in the study 
area to those in the control area indicates that total 
production in the area of the sewage field was approx- 
imately double that in the control area. Although 
enrichment by the sewage undoubtedly contributed to 
this increase, it may be speculative to assign all of it 
to that cause. The submarine canyon off Newport, with 
its upwelling of nutrients from aphotic depths, may 
well have contributed to the difference. 

The relationship of red tide to waste discharge is 
still doubtful. Its presence in the study area and its 
general absence in the control area may reflect chance 
conditions encountered in the relatively few cruises to 
the control area. Alternatively, the net water movement 
through the study area requires less than 24 hours and 
the previous history of the water is of some importance 
in determining its phytoplankton population. In addi- 
tion, red tide was commonly heavy upstream and 
offshore from the sewage field. This, together with the 
fact that heavy red tide occurs in coastal areas quite 
remote from any outfall, makes it doubtful that sewage 


TABLE 9 


The mean daily production at the surface for the four seasons in 
the study area off Orange County, and the control area off Camp 
Pendleton. The data are expressed as metric tons of carbon fixed 
per day over a total area of 60.2 square kilometers or 17.6 square 
miles. Figures in parentheses are from samples collected in the study 
area during periods of red tide. 


Season Study area Control area 
Winterloe eae ee 4.8 (7.9) No observations 
Spring 22. cee 9.7 (11.9) 4-2 
Summers > .ee 2 10.0 (no red tide) No observations 
Kall.c2 scsus see ee 4.4 (23.6) 2.4 


is related to the periodic development of heavy red 
tide along this coast. 


4. FACTORS POSSIBLY RELATED 
TO OBSERVED CHANGES 


From the work in fresh water, Hynes (47) made 
certain assumptions regarding the sequence of changes 
that may be of biological significance in receiving 
waters; a) dissolved oxygen content will be sharply 
reduced at the point of sewage discharge and will rise 
gradually downstream as the BOD is satisfied; b) the 
addition of suspended and dissolved matter in the 
effluent will sharply decrease water transparency near 
the outfall and return to normal downstream; c) the 
suspended material will settle out forming a graded bed 
of sediment along the course of the current; d) toxic 
and inhibitory substances normally absent, or present 
in low concentrations, together with e) nutrient sub- 
stances normally present in small quantity, will be in 
high concentration in the vicinity of the outfall and 
will gradually decrease downstream as both chemical 
and biological processes remove, dilute, or alter them. 
Each of these factors as they affect productivity in 
open coastal waters are discussed as follows: 


a. Dissolved oxygen— 


The oxygen content of the sewage field rapidly ap- 
proaches that of the background waters, with the most 
rapid increase occurring during initial mixing when the 
plume rises to the sea surface. At the surface, oxygen 
is present at about 75% of saturation. Subsequent eddy 
diffusion brings about a further increase. This rapid 
mixing rate with receiving water which contains oxygen 
at 90-140% saturation more than satisfies BOD re- 
quirements and the oxygen pattern within the spread- 
ing field is entirely different than that described by 
Hynes (loc. cit.). It is unlikely that levels of oxygen are 
sufficiently low to interfere with the respiration of 
pelagic organisms within the surface field. 


b. Transparency— 


Low transparency of the water near the boil obvi- 
ously limits the euphotic layer to an extremely shallow 
depth. Over the area as a whole, turbidity is of great 
significance to the benthic biota but is of less significance 
to the pelagic life of the area. Due to the variability in 
direction and magnitude of the currents, unlike the 
confined flow in streams and rivers, steady state condi- 
tions do not exist at any point in space except in the 
immediate vicinity of the boil. Consequently the portion 
of the total area that is continuously affected is rela- 
tively small. Dilution at a relatively high rate also 
restricts the magnitude of the effect. 

In considering pelagic populations, the situation is 
comparable to that in streams and rivers if one deals 
with events along the trajectories of the sewage field, 
except for the important difference that the availability 
of the diluting water and its contained population is 
very large and the rate of dilution is high and con- 
tinuous. 

Direct measurements of water transparency indicate 
that at and near the boil the transparency of the surface 
water was reduced to 10% or less of the transparency 
in air (Kahl transparency meter, path length = % 
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meter). At a depth of about 5 to 10 meters transparency 
was very variable depending on whether or not the 
meter was in the rising sewage plume. At about 2 hours 
downstream, transparency at the surface was about 
70-75%, and at 5 meters or more it was about 80-85%, 
a value which was approximately that of background. 
Further downstream the surface transparency increased 
slightly but somewhat lower transparency was encount- 
ered beneath the surface. The averaged data are pre- 
sented in Table 10. Inspection of these data indicates 
that alteration of natural conditions is relatively short- 
lived. The sinking of the band of lowered transparency 
to at least 10 meters after 6 hours, and its increasing 
vertical dimension, seem to reflect vertical diffusion 
combined with settling of particulate material. 

The effect of this reduction in water transparency on 
photosynthesis of phytoplankton may be estimated by 
the ratio of production to productivity of subsurface 
samples. Such a comparison is made in Figure 51. No 
red tide was present during these cruises. Samples were 
collected at 2-hour intervals in the path of a dye marker 
or drogue released in the boil. Obviously a blanket of 
red tide would produce the same effect as a loss of 
water transparency due to sewage but the data would 
not be directly attributable to the discharge of wastes 
into the waters. Inspection of this Figure shows that at 
8-12 meters the effect of light reduction is apparently 
negligible after 12 hours. There is no significant effect 
on productivity apparent at this depth that could be 
attributed to the presence of sewage. At 1-3 meters 
the production and productivity relationship is dis- 
turbed by the presence of sewage during the first 10 
hours. Beyond that time the level is similar to that 
from the deeper samples. Thus, the effect of reduced 
light transparency on phytoplankton production is 
relatively short lived. This agrees with Hynes’ conclu- 
sion, although the time and distance scale will vary 
depending on the particular environmental situation. 


c. Sedimentation— 


The formation of graded beds of debris deposited by 
sedimentation of particulate matter from the sewage 
is of minor importance to the pelagic crganisms with 
which this project is concerned. Therefore no studies 
were made of rates of settling or deposition. It is worth 
noting, however, that although the bottom may be 
greatly altered ecologically (Hartman (44) ) there is 
no observational evidence that static beds of material 
continuously accumulate despite many years of opera- 
tion of various marine outfalls along this coast. 


d. Toxic and inhibitory substances, and 
e. Nutrients— 


These two headings are dealt with jointly since the 
methods employed do not determine the specific physio- 
logical role of the various substances for which analyses 
were made. For example, many nutrients can inhibit 
photosynthesis if they are present in too high a con- 
centration. Our data show that after an initial depres- 
sion, productivity starts to rise in the sewage field. 
Continued dilution results in a return to background 
levels of all nutrients and other parameters which were 
measured. Other processes which occur at the same 
time as physical dilution can result not only in a further 


TABLE 10 


Average transparency of water at stations occupied at timed 
intervals along the course of a diluting sewage field. The data are in 
percent transmittance along a light path of one-half meter compared 
to the transparency in air. 














Hours 
0 2 4 6 
Depth (Boil) 
Surfacessoee =e. Bf 9% 75% 79% 81% 
5-10 Meters________ 66% 84% 70% 70% 
10+ Meters___.____ 62% 82% 75% 61% 





alteration of absolute concentrations of substances but 
also in a change in their relative composition. These 
changes may be biological (uptake and release of sub- 
stances by living organisms such as bacteria and 
phytoplankton), physical (adsorption on particles), and 
chemical (breakdown and synthesis of compounds). 
Thus, the pattern of productivity found in a diluting 
sewage field can be hypothetically attributed to the 
change in concentration, or to the alteration in com- 
position of the effluent, or to a combination of the two. 
The analytical methods in use during the course of 
this study have been aimed at determining the con- 
centration of selected compounds generally considered 
to be of biological significance. 

The hypothesis that the sequence of events described 
above for productivity can be related to changing 
concentration of sewage was tested by the experimental 
exposure of a natural population to various concentra- 
tions of filtered effluent. The results of this experiment 
are listed in Table 8 of the Annual Report for 1962-63 
(1) and indicate that high concentrations of effluent as 
found in the boil were indeed inhibitory. More dilute 
effluent stimulated productivity. It is worth noting that 
these effects were still evident after 14 hours storage of 
the mixture of phytoplankton and effluent. 

Attempts to relate biological uptake of nutrient salts 
to reduction in concentration of the same nutrients in 
the medium have not been conclusive. In one field 
experiment, changes in the concentration of nitrate, 
phosphate and silicate were measured in both the light 
and dark bottles used for productivity determinations. 
The changes which occurred between the beginning and 
the end of the incubation period were less than the 
limits of error of the analytical procedures. 

The three major nutrient substances for which ex- 
tensive analyses have been made are nitrate, phosphate 
and silicate. These three, because of their biological 
importance, have been among the most widely mea- 
sured and considered. The first two, nitrate and phos- 
phate, are universal requirements for algal growth, and 
the last, silicate, is essential to the diatoms. The corre- 
lation between productivity and these nutrient salts 
in the field was poor whether the data were considered 
as a whole or for each individual cruise. Figures 52 
through 60 illustrate the problem and show both the 
productivity and chemical data for different cruises 
plotted as a function of time of drift from the boil. The 
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concentration of nutrients was maximum at the boil 
and fell off with distance or time following a relatively 
simple curve. Productivity, whether considered as a 
function of water volume or of units of chlorophyll, 
was low near the boil, rose to a maximum value and 
then fell off. Although the basic pattern was the same 
for all cruises, the slopes of the lines and the relative 
values are different for different cruises. 

Since the basic pattern remained the same for the 
different cruises an arbitrary system of zones was set 
up based on the productivity per unit of chlorophyll. 
These zones, 5 in number, can be characterized as 
follows: 1) the boil and immediately around it where 
the productivity was essentially zero; 2) the area down- 
stream from the boil in which the productivity per unit 
of chlorophyll was low; 3) the area downstream from 
zone 2 in which the productivity per unit of chlorophyll 
started to recover from the initial inhibition associated 
with sewage; 4) the area downstream from zone 3 in 
which productivity per unit of chlorophyll reached 
maximum values; 5) the area downstream from zone 4 
in which productivity per unit of chlorophyll decreased 
toward background levels. 


TABLE 11 


Correlations between surface productivity per unit of chlorophyll 
and nitrate, phosphate and silicate in various zones around the 
Orange County Treatment Plant outfall. The numbers in parentheses 
refer to the number of samples. 





Zone 








Productivity and Nitrate 


10-62 Sere ees Sete Be ee .16(12) 

11—-6225_ - .71(3) 

12-022 See |e eee ee oe ee .43(12) 
9-63 2 he. — .09(5) .12(5) .61(8) 

11-6322, o> — .52(3) — .52(31) 

12-60 Se ee gre ore [ee ee tee .32(6) .90(5) 
S04 ee | a ee eee SPAM Gn? — .98(5) fF .85(6)* 
O-0452 22 eee eee ae | ee — .60(19) TF 


10-02 221223 | See ee See er ete — .02(12) 
Mao 22 See £94(3)\ fe eee we — .00(6) 
12-02 228 2 ees eee ee .65(12)* 
3-0325 ste ae eee .98(5) TF .59(3) 
4-63.52 sees Se ee eee — .16(9) .00(4) 
URES Ree a .02(3) 
12-63.) sie ee eee ee eee — .32(7) .37(5) 
nig CM | es, Sa apes .15(7) .75(5) .41(6) 
G-642 2 32 |te See ees |S eae eee .01(21) 


10-62 (Fea eee ee ee oe — .05(12) 

t—b2ee oe $48(3) |S e a ee — .64(6) 

1eH620) 22 Soa ee oo |e eee — .51(12) 
3—Oe =e —=§.98(5) fF — .97(6) Tf — .50(8) 

11-632 — 25 0.00(3) 

12-63... 3-|2 028 LS Se — .53(7) .11(5) 
O-O4) 3. 221k ee ee — .44(7) — .76(5) . 25(6) 
6-64... |e See coe ee .14(21) 


* Significant 
+ Highly significant 


The statistical relationship in each zone between 
productivity per unit of chlorophyll and each of the 
nutrients was computed separately for each cruise. 
Data believed to be of questionable validity for tech- 
nical reasons (instrumental failure, etc.) were not 
considered. Treating the data in this fashion and using 
the data from each cruise separately, reduced the 
sample numbers for each calculation but improved the 
chance of detecting significant relationships. Correla- 
tions for the data from zone 1 were not calculated since 
the productivity in this zone was consistently around 
zero, although the nutrient concentrations were high 
and variable. The consistently low levels of production 
in this zone suggest either that the levels of nutrient 
present were high enough to be inhibitory, or specifically 
inhibiting factors were present. The data do not indicate 
the threshold level of concentrations at which inhibition 
would take place. In the other zones, the correlations 
found for the data from various cruises are listed in 
Table 11. These data indicate that from cruise to cruise 
and from zone to zone there was tremendous variation, 
although high or significant values occasionally appear. 
The lack of general statistical relationships indicates 
that these nutrients in the area are not limiting. 

Similarly treated data from the control area are 
presented in Table 12. The zones, in this case, corre- 
spond to bands which roughly parallel the shore and 
which are numbered consecutively from the shore 
outward. This arrangement takes into account the 
expected decrease in production with distance from the 
shore (Doty & Oguri (67) ). The correlations are ex- 
tremely variable and, in general, not high except in 
zone 4, the furthest offshore. This zone was occupied on 
only one cruise. Apparently the nutrients determined 
here do not normally limit production in inshore waters. 


TABLE 12 


Correlation between surface productivity per unit of chlorophyll 
and nitrate, phosphate and silicate in the control area near Camp 
Pendleton. The zones in this case are numbered consecutively from 
inshore to offshore. The numbers in parentheses refer to the number 
of samples. 











Zone 
Date 1 | 2 3 4 
Productivity and Nitrate 
10-62____- 0.95(3)) fee wee ee 0.68(4) 
12-6362 502) eee oe ee .33(3) 
3-64_____ — .44(11) .08(13) — .36(4) — .95(3) 
GH642 22.2 eation © Cee eee eee oe .52(4) 
Productivity and Phosphate 
10-62____- Ei 2(3) ew ee ener — .78(4) 
12-6353. | A .94(3) — .35(6) 
3-64____- .02(11) — .06(13) .42(4) — .86(3) 
HAM RPT! bol Pt eke | te oe ee — .22(4) 
Productivity and Silicate 
10-62____- 299(3)* [peewee — .08(4) 
NW Oe he en CN Ri. ce .70(3) .72(6) 
3-64_____ .90(11)F .35(13) . 25(4) — .90(3) 
G64 ei eae ae .87(3) — .02(4) 


* Significant 
+ Highly significant 


D. ULTRAVIOLET ABSORPTION AS A MEANS OF TRACKING 
A SEWAGE FIELD 


One of the problems arising out of our work in the 
marine environment is that of physically tracing the 
course of the currents carrying the sewage and of as- 
sessing the relative concentration of the sewage at the 
time and place of our stations. The use of current 
markers, such as drogues or dye patches released at the 
boil, can effectively give information on the course 
followed by the sewage. However, at any given point 
along the track of the current marker there is no really 
convenient way of ascertaining that the stations occu- 
pied were in the most concentrated portion of the more 
extensive effluent field. 

The obvious solution, dyeing the entire effluent field, 
was tried once. Although successful, the quantity of 
dye required made this method too expensive for 
routine use. 

Analyses for components of the sewage not normally 
present in sea water or present in very small quantity 
can yield the desired information, although some are 
not primarily used for that purpose (BOD and coliform 
counts). If the effluent is of low salinity, salinity de- 
terminations can also be used. These analyses, however 
are neither sufficiently sensitive nor rapid enough for 
field use. Analysis for ammonia, normally present in 
sea water only in very low concentrations, has been 
found to be useful for some purposes although it is 
inconvenient as a field method on board ship and it is 
not adequate at high dilutions of sewage. 

In previous reports we have described the finding 
that high optical absorbance of Millipore-filtered water 
to light at a wavelength of 230 millimicrons is related 
to the concentration of sewage. This method is fast 
and simple, it requires less than 100 ml. of filtered water 
from the sample, and direct reading of the absorbance 
can be made in a spectrophotometer. The Beckman 
Model DB has been found to be an ideal instrument 
for this purpose since it is a double-beam instrument 
and the background reading is automatically subtracted. 

The results of routine measurement made during the 
last half of this project are plotted in Figures 61, 62, 
63, & 64 for the different seasons and in Figure 65 for 
all seasons together. The distribution of absorbance 
shows the expected high values around the boil and a 
decrease toward background values with increasing 
distance from the boil. Differences in shape of the 
affected areas for the different seasons reflect the differ- 
ences in currents expected for the different seasons. 

The pattern of high absorbance near the boil and a 
decrease with distance was similar to that found for 
the chemical substances, notably ammonia. Since am- 
monia may also indicate the presence of sewage, this 
suggested that the relationship between the two could 
serve as a preliminary measure of the adequacy of the 
UV method. Table 13 lists the correlations and regres- 
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TABLE 13 


Correlation and regression between ammonia and ultra violet 
absorbance, Orange County outfall, California. 





Number of 
Date samples Correlation Regression 

9-11-63________- 12 0.524 114.8 
11- 9-63________- 26 0.475 106.7 
11—-10-63_________ 18 0.924 394.2 
12-14-63_________ 25 0.760 187.5 
2— 8-64___...._- 34 0.853 170.0 
6-17-64________-_ 28 0.974 231.1 
Totaiewt #14. 143 0.732 263.4 





sions found during each of six of our cruises for which 
comparable data are available. The totaled data for 
all cruises are also included. Correlations were high in 
all but the earliest cruise, and the variation in regression 
from cruise to cruise is evident. The high correlations 
indicate that both UV and ammonia may be functions 
of the same process, i.e., dilution, and the regressions 
indicate that probably they are not direct functions of 
one another. 

Data for cruises in the sewage-free control area off 
Camp Pendleton for which comparable information 
exists are presented in Table 14. There is a significant 
correlation in only one run of 4 listed. Even though 
UV absorbance and ammonia were so low that the 
limits of sensitivity of the methods may have contrib- 
uted to the observed variations, comparison with data 
from the study area shown in Table 13 indicates that 
the UV absorbance in the area where sewage is absent 
is not correlated with ammonia. 

Since both elements (UV absorbance and ammonia) 
show good statistical relationship to one another within 
a sewage field but not outside and presumably are 
functions of the concentration of sewage, it would be 


TABLE 14 


Correlation and regression between ammonia and ultra violet 
absorbance data found during cruises to the control area off Camp 
Pendleton, California. 














Number of 
Date samples Correlation Regression 
5-63.22 ee 7 0.090 8.4 
12-6350 eee 15 0.700 29.3 
3-645) a Soe 5 0.662 2.3 
6-04 50574 fees 5 — 0.357 23.6 
Total s30ees- 32 0.286 34.1 
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reasonable to expect similarly good correlation between 
them and the degree of dilution of the sewage field. The 
correlations in Table 15 show that the expected good 
relationship was not found. The diffusion data for these 
correlations were obtained from the overall diffusion 
equations for each particular run and the maximum dye 
concentration (Cmax) was computed for the time or 
distance at which the stations for UV and ammonia 
were occupied. This was necessary since the field ob- 
servations for diffusion data are collected independently 
and by a different method from the other information. 
The diffusion data represent the overall behavior of a 
relatively large dye patch in a larger sewage field. The 
chemical and biological data represent the findings at 
a series of discrete locations just ahead of the dye patch. 
Other differences in methods of treating the data also 
may be among the reasons for the apparent lack of 
good correlation. One limitation may be associated with 
the attempt to use dilution equations for computing 
the instantaneous concentration at specific points of 
field observation. It is also possible that UV absorbance 
and ammonia are not conservative properties. Testing 
on these points will continue. 


TABLE 15 


Correlations of ammonia and ultra violet absorbance with calcu- 
lated maximum concentration of rhodamine-B dye (Cmax) at the 
time or distance from the boil corresponding to sampling stations for 
the chemical samples. The field observations were all made in the 
study area off Orange County, California. 


Cmax—distance Cmax—time 

Date UV NHsz UV NH; 

9-63 24: Le aes. 0.333 0.958 0.725 _— 
Tl-63 2.0". Se 0.968 0.697 — 0.968 — 0.858 
n= wl ES 8: jn a: n=83 
6-64. -2u esa fabs 0.554 — 0.228 0.449 0.261 
n= 4 n=4 n=4 n=4 
6-642 2 ane 0.408 0.801 0.019 0.804 
n= n= 6 n= 6 n = 6 
Totals ee 0.634 0.390 — 0.028 — 0.193 


n= 13 
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E. CHEMICAL PROCEDURES 


1. METHODS 


Sampling: Water samples were collected from the 
surface with the aid of buckets and from subsurface 
with the aid of a snatch bottle, Nansen bottles, Van 
Dorn type sampler, or a new subsurface sampler de- 
veloped during the course of this project (Fay and Smith 
(68) ). All samples for biological and chemical studies 
were collected with inert plastic samplers or plastic 
buckets. Subsurface samplers were lowered with the 
aid of a bethythermograph winch to depths of 200 m. 
or less and with a hydrographic winch to depths below 
200 m. 

Filtration: Subsamples requiring filtration were passed 
through a Millipore filter, 47 mm. in diameter, of a 
nominal porosity of 0.8 microns. Rapid filtration was 
accomplished with the aid of a Millipore vacuum- 
pressure pump. 

Temperature: Determined to +0.1°C. with the aid of 
a bucket thermometer, a thermometer mounted inside 
of the plastic subsurface samplers, a bathythermograph, 
or to +0.02°C. with reversing thermometers. 

Salinity: Determined by titration with silver nitrate 
to the chromate end point (Strickland and Parsons 
(69) ) or to the fluorescein end point (Miyake (70) ) to 
+0.02°/oo or with an inductive conductivity meter to 
+0.003°/ 9. Samples were stored in polyethylene bottles 
in a dark, cool location and analyzed within two weeks 
to one month after collection. 

Dissolved Oxygen: Determined by the Winkler titra- 
tion as modified by Pomeroy and Kirschman (71) with 
a precision of +0.05 ml./L. Determination of the per- 
cent saturation of dissolved oxygen was based upon 
the predicted concentrations of dissolved oxygen in sea 
water as computed by Whipple and Whipple (72). 
Thiosulfate solutions were standardized with potassium 
dichromate as recommended in H.O. Publication No. 
607 (73). Samples for the determination of dissolved 
oxygen were taken directly from the subsurface sam- 
plers or from the bucket containing a sample of surface 
water immediately following collection. 

pH and Alkalinity: A subsample of 100 ml. of unfil- 
tered sea water was aliquoted with the aid of a gradu- 
ated cylinder. The pH of this subsample was determined 
with the aid of a Beckman Expanded Scale pH meter 
to +0.01 pH units. Hundredth normal hydrochloric 
acid was added to the sample whilst stirring until a pH 
of 4.00 was obtained. The milliequivalents of acid added 
to bring the pH of the sample to 4.00 was taken as an 
index of the alkalinity of the sample to +0.02 milli- 
equivalents/L. 

The following analyses were routinely performed 
assuming the sample to contain only ‘soluble’ com- 
ponents (those passing through a Millipore filter of a 
nominal porosity of 0.8 microns). The 400-600 ml. of 
filtrate thus obtained was subaliquoted in volumes of 
50 or 100 ml. into 4 oz. polyethylene bottles with the 
aid of graduated cylinders of these respective capacities, 
frozen, and maintained frozen at about —8°C. until 
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thawed and brought to room temperature for the de- 
termination of the particular analyte of interest. Re- 
agents were added directly to the plastic bottle contain- 
ing the subsample with the aid of Miscomatic Dispensers 
and the reactions resulting in the formation of the 
colored complexes were allowed to come to completion 
in this container at room temperature or in a water 
bath (nitrate only). Occasionally the analyses were 
performed upon filtered samples aboard ship immedi- 
ately after collection with essentially the same results 
obtained with fresh as with frozen samples (similar 
studies have also been reported by Stefansson and 
Richards (74) and within the expected limits of pre- 
cision of the analysis. 

Phosphate, ‘‘Reactive’: Determined by the methods 
described by Strickland and Parsons (69). Replicate 
50 ml. samples agree with a precision of +2% and with 
an accuracy of +2%. 

Silicate, “‘Reactive’’: Determined by the method de- 
scribed by Chow and Robinson (75) and defined by 
Strickland and Parsons (69). Results are reported with 
a precision of +5% and an accuracy of +10%, for a 
sample of 50 ml. 

Nitrite and Nitrate: Determined by the methods de- 
scribed by Strickland and Parsons (69) with the follow- 
ing modifications. A sample of 100 ml. of filtered water 
was transferred to an opaque, 4 oz. polyethylene bottle 
and the appropriate reagents added directly to the 
sample in this container. Precision varies with the con- 
centration of the analyte being +10% in the range of 
0.2 to 2 microgram-atoms of NO;3-N/L. and +5% in 
the range from 2 to 40 microgram-atoms of NO3-N/L. 
and +5% in the range from 2 to 40 microgram-atoms 
of NO;-N/L. Nitrate was not corrected for nitrite as 
nitrite was not observed to exceed 5% of the nitrate 
found in any one sample in which both were deter- 
mined. No difficulties in the analysis were experienced 
as a result of the substitution of plastic bottles for 
glass flasks. 

Ammonia: Determined by Nesslerization with the 
following modifications. All reagents were made up in 
double distilled, deionized water. Sodium tartarate 
replaced sodium, potassium tartarate. Sodium tartarate 
was recrystalized from 0.05% sodium hydroxide by the 
addition of one volume of acetone to five volumes of a 
saturated, alkaline solution of sodium tartarate. Ness- 
ler’s reagent was prepared according to the directions 
provided in ‘Standard Methods for the Analysis of 
Water, Sewage, and Industrial Wastes’’, 10th Ed. (75). 
Each sample of 50 ml. of filtered sea water received 2 
ml. of 30% tartarate, 3 ml. of 20% sodium hydroxide, 
and 1 ml. of Nessler’s reagent. The timing of the addi- 
tion was closely controlled. The addition of the sodium 
hydroxide was followed within 1-5 minutes by the 
addition of the Nessler reagent. The color was allowed 
to develop for 7-10 minutes and the optical density 
was determined at 420 millimicrons in a cuvette with 
a 10 cm. light path. Precision varies with the concen- 
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tration of ammonia being +20% in the range of 0-2 
microgram-atoms of NH;-N/L. and +5% in the range 
of 2-100 microgram-atoms NH;-N/L. 

General Remarks: In all analytical procedures which 
are subject to a correction factor related to the salinity 
of the sample, standard solutions were made up in local 
sea water deficient or depleted in the specific analyte 
(e.g. ammonia, nitrate) or with artificial sea water of 
a composition (phosphate, silicate) similar to that of 
local sea water. Recovery samples were also included 
with the standards by adding known amounts of the 
specific analyte to duplicate samples. 

Polyethylene bottles were color coded for the various 
components determined during this study by spray 
painting the sample bottles. As each bottle was filled, 
it was labelled with the appropriate station number, 
the depth from which the particular sample was col- 
lected, then frozen immediately in the ship’s deep freeze 
and kept frozen until analyzed in the shore base la- 
boratory. 


2. COMMENTS ON THE ANALYTICAL 
PROGRAM 


In keeping with the objectives of the demonstration 
grant to develop and evaluate techniques for the study 
of the biological effects of waste disposal in the marine 
environment, the analytical program is reviewed with 
respect to the potential value of the methods and 
techniques used. 

First, the analytical program was primarily designed 
to support the dilution and productivity studies, but 
some information was also obtained on the variation in 
chemical composition of an effluent field. Diffusion 
studies required knowledge of the vertical stability of 
the water column both inside and outside of a field of 
effluent waste as vertical stability is related to both the 
vertical and horizontal diffusion of the waste field. 
Independent determination of the patterns and rates 
of waste dilution also was required to verify the con- 
clusions obtained from the studies of the dye diffusion. 
This was accomplished by quantitative analyses of 
potentially useful chemical tags. Primary planktonic 
productivity varied markedly with the composition of 
the medium, and an attempt was made to relate plank- 
tonic productivity with variations in the chemical 
composition of both the diluting effluent field and the 
receiving waters. Also, differences in composition were 
established in order to distinguish the sewage field from 
background receiving waters. Finally, the receiving 
waters themselves deserved attention in order to 
establish normal background conditions as compared 
to the sewage field resulting from essentially steady- 
state discharge of effluent into these waters. 

The significance of each of the determinations per- 
formed on a routine basis to each of the portions of 
this study shall now be reviewed as a guide to future 
work of a related nature. Figure 66 is a diagram of the 
sampling scheme. At some stations, not all of these 
determinations were carried out, but a selection was 
made according to specific requirements. 

Temperature is the most easily determined of all the 
variables measured in this study. Knowledge of the 
distribution of temperature vertically through the 
water column was required to determine the density 


structure and stability of the water column. Data on 
vertical stability was important in studies of vertical 
coefficients of eddy diffusion and of the vertical dis- 
placement of horizontal layers of the water column by 
internal waves and other physical processes. Tempera- 
ture also may be an indicator of the horizontal distribu- 
tion of water masses and of the solubility of various 
substances, particularly gases, in sea water. It is often 
related to the distribution patterns of plankton, nekton, 
and fixed benthic organisms. Temperature is not 
especially useful in defining a field of waste other than 
the unique case of following a discharge of hot water 
from an industrial plant. 

The horizontal and vertical distribution of salinity 
may be used as an indicator of the distribution of water 
masses and of the associated distribution of various 
organisms in the sea. It is also related to the activities 
of solutes in sea water. Salinity may be determined 
relatively rapidly to a high degree of precision with the 
aid of a conductivity meter. Caritt (76) has recently 
reviewed the application of the conductivity meter to 
this determination and has pointed out the advantages 
of this instrument. Salinity is as useful at the level of 
sensitivity of the determination, which for practical 
purposes may be taken as +0.005°/o, as are any of 
the other indicators of the presence of effluent wastes 
routinely determined in this study (phosphate, silicate, 
and ammonia). If sufficient data are taken to establish 
the limits of variability of the receiving waters, salinity 
may be the most convenient determination available 
with which to estimate the initial dilution of sewage 
and the short-term pattern of subsequent diffusion. 

Dissolved oxygen content is an important constituent 
in assessing the quality of the receiving water and the 
effects resulting from effluent discharge (see Reference 
(75) ). Oxygen has been found to vary in concentration 
in a significant, seasonal pattern in the nearshore waters 
of southern California to a depth of about 200 meters. 
Below 200 meters, in the deep basins of the southern 
California continental borderland, dissolved oxygen 
decreases to a minimum of 0.2 ml/L (Emery, (77) ). 
Oxygen does not serve as a distinctive tag of the effluent 
field at the Orange County outfall for the following 
reasons: 1) dilution of the wastes proceeds so rapidly 
that the concentration of dissolved oxygen approaches 
saturation within one-half mile of the boil, 2) dissolved 
oxygen has not been present at less than 75% of satura- 
tion anywhere in the field of effluent waste with the 
exception of the boil where it is equal to or greater than 
70% of saturation, 3) in an effluent field floating on 
the surface of the receiving waters, equilibration of 
dissolved oxygen with oxygen from the air may be 
expected to rapidly restore any localized deficiencies of 
this gas, and 4) submerged fields established over out- 
falls with efficient diffuser ports (as at the Hyperion 
outfall in Santa Monica Bay) are characterized by a 
high initial dilution of 1:100 or more, and differences 
in oxygen content between the field and background 
waters become insignificant. 

For the purposes of this study, the determination of 
dissolved oxygen by the Winkler method has proved 
to be completely adequate. Titration is simple, reliable, 
essentially fool-proof, and free of the problems en¢ount- 
ered with oxygen electrodes now available. Two different 
oxygen electrodes were obtained during the course of 
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this program from commercial sources. Neither func- 
tioned adequately for the purposes of this study. Until 
a reliable, stable, 7m situ oxygen electrode capable of 
submersion to 200 m. is developed, the Winkler titration 
will continue to be the more dependable and useful 
technique for determination of dissolved oxygen in sea 
water. Winkler reagents, together with appropriate 
titration apparatus, should be available in reserve for 
any electronic equipment when field surveys are con- 
ducted. 


Use of the Beckman Expanded Scale pH meter aboard 
ship permitted rapid and reliable determination of the 
pH of unfiltered samples of sea water. The pH decreased 
with depth and varied seasonally. Titration with 0.100 
N.HCIl to determine the alkalinity of the sample of sea 
water is accomplished with about the same precision 
as the technique employed to obtain the sample aliquot 
from which this variable is determined (+1%). The 
alkalinity of surface samples increased about 20% 
above background in the immediate vicinity of the boil. 
Alkalinity decreased slightly with depth, and varied 
with season and with water mass. Determination of 
alkalinity was not useful in following a diluting field 
of effuent waste because the precision of the determina- 
tion was less than the differences observed between 
background and the sewage field at radii greater than 
(about 1 mile) from the boil at the Orange County 
outfall. Alkalinity data can be useful in following water 
masses and in determining the productivity of surface 
waters on the basis of changes in pH, carbonate con- 
tent, and dissolved oxygen (Sargent and Hindman, 
1943). 


“Reactive” phosphate is one of the most interesting 
components quantitatively determined in this study. 
It is reliably determined although the ease and precision 
of the analysis could be improved by using reducing 
agents and molybdate which are specifically purified 
to remove phosphate, and by developing reducing 
agents which are stable for at least 30 minutes. Phos- 
phate serves as a reliable tag for the presence of effluent 
wastes provided a sufficient number of background 
samples are collected at the same time to provide a 
clear distinction between the effluent field of waste and 
the receiving waters. Phosphate also serves as a tag 
for the study of the dynamics of inshore waters and the 
identification of different water masses. Phosphate is 
traditionally determined as an indicator of the fertility 
of the waters in the euphotic zone. 


“Reactive” silicate likewise serves to identify the 
extent and dilution of an effluent field if sufficient 
synoptic determinations are made to establish the im- 
mediate composition of the receiving waters. Silicate 
can be determined with sufficient ease and precision to 
be a useful indicator of the initial effluent field and of 
the distribution of the water masses found in the in- 
shore area. Silicate appears to be unsatisfactory as an 
index of the fertility of the receiving waters of Southern 
California. 

At this time, ammonia is the most useful chemical 
tag of those with which to describe the extent and 
dilution of an effluent field. As useful as ammonia has 
proved to be, this substance is also the most elusive 
analyte with respect to the ease, accuracy, and de- 
pendability of the procedure. A rapid, reliable colori- 


metric method for the determination of ammonia in 
sea water is still to be developed. Nesslerization offers 
the promise of such a method but, for as yet unknown 
reasons, it is not always dependable when applied to 
samples of sea water. Since a very large number of 
samples are necessary adequately to describe the dilu- 
tion and distribution of a field of effluent waste, either 
colorimetric or potentiometric techniques will be re- 
quired to expedite this analysis. 

Nitrite is determined with comparative ease and 
accuracy; however, with the possible exception of a 
detailed study into the rate of oxidation of ammonia 
through the intermediate stage of nitrite to nitrate, 
nitrite offers little advantage as a tag for the presence 
of the field of effluent wastes or as an indicator of the 
dynamic processes of the inshore water masses. 


The nitrate analysis is an inconvenient one; when 
nitrate is determined by the method described by 
Stickland and Parsons (69), a reduction time of about 
20 hours must be sustained before color development 
is completed and the method therefore is ill suited to 
routine use in the field. Other methods, such as the 
polarographic reduction described by Chow and Robin- 
son (75) or the procedure for the determination of total 
dissolved inorganic nitrogen described by Riley and 
Sinhaseni (79), while not satisfactory for use aboard 
ship, may provide useful cross checks upon the accuracy 
of the other individual determinations by colorimetric 
techniques for nitrate, nitrite and ammonia. 

For accurate determination of dissolved substances, 
filtration is indispensable when processing samples from 
turbid inshore waters. A filter with a nominal porosity 
of 0.8 microns was used in this study because the amount 
of suspended material found in the inshore waters was 
sufficiently great to rapidly plug filters of a finer porosity 
when filtering volumes of 0.5 to 1.0 liters. 

Rapid sampling through the water column is essential 
to the study of natural waters. An improved sampler 
was developed which has increased the rate of sampling 
about 50% over that formerly obtained with the 
Van Dorn type of sampler, and is a more rugged, 
seaworthy device. 

More reliable and rapid analytical procedures are 
needed for ammonia, nitrate, total fixed inorganic 
nitrogen, phosphate, and various dissolved organic 
compounds found in sea water and sewage. The Beck- 
man DU spectrophotometer should be replaced with a 
more seaworthy instrument of comparable or increased 
sensitivity and stability. A simple split beam filter 
photometer of appropriate sensitivity could well solve 
the operational shortcomings of the DU spectropho- 
tometer (sensitivity to humidity, clumsy cell compart- 
ments, restricted sample cells, bulk, etc.) at a con- 
siderably lower cost in instrumentation. Standardization 
of analytical procedures is also needed. Precedent has 
been established in this direction by the formulation of 
‘“‘Recommended Interim Procedures for Measurements 
in Biological Oceanography” (80). 

The most useful properties or components for de- 
scribing background processes and defining the sewage 
field have included temperature, salinity, dissolved 
oxygen, pH, phosphate, silicate, ammonia and nitrate. 
Quantitative determinations of alkalinity and nitrite 
were of less value although they may be of specific 
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interest for considerations not included in this study. 
The principal contribution of the chemical program lay 
in describing the distribution of a relatively large 
number of variables on a routine basis throughout a 
two-year period, which were related to studies of 
primary planktonic productivity, to the dilution of a 
field of effluent wastes, and to the physical and biolog- 
ical properties of the receiving waters. 


3. SOME ASPECTS OF THE CHEMICAL 
COMPOSITION OF BACKGROUND 
WATERS 


In order to support the productivity and dye diffusion 
studies, a flexible pattern of chemical sampling was 
established within a two mile radius of the Orange 
County outfall to permit an examination of the com- 
position of the receiving waters as well as that of the 
diluting sewage field. Since the direction of flow of the 
sewage field varied, the location of the stations describ- 
ing background composition also changed. Primary 
indicators useful in detecting pollution were found to 
be silicate, phosphate, ammonia, and salinity. Ancillary 
information was also obtained using pH, dissolved 
oxygen, and nitrate data. Generally, the studies of the 
effect of sewage wastes upon planktonic productivity 
concentrated attention in the sewage field; thus there 
were 289 stations in the sewage field and 141 which 
were taken to describe the background concentration 
of the receiving waters. In addition, a control station 
was routinely occupied in about 100 fathoms of water 
(98 hydrocasts during two years) over the axis of the 
Newport Submarine Canyon about 2 miles southeast 
of the boil. Other control areas at various locations 
along the southern California coast and near Catalina 
Island some 20 miles offshore were examined occasion- 
ally for productivity, dye diffusion, and chemistry. 

Background samples were collected through the 
water column generally at the standard depths of 0, 
5, 10, and 25 meters or in conjunction with the studies 
of planktonic productivity, at 0, 3, 6, and 12 meters. 
This pattern of sampling provided adequate data with 
which to outline the instantaneous status of the vertical 
composition of the receiving waters; repetitive sampling 
through the year revealed successive changes in com- 
position at selected depths through the water column. 

Background composition of the receiving waters near 
the Orange County outfall from mid-1962 to mid-1964 
is depicted in Figures 67 and 68. Seasonal trends are 
evident in spite of certain irregularities in the curves. 
Temperature showed the expected minima in the late 
winter and spring and maxima in the summer months. 
The salinities were also high in the summer but were 
generally lower in the early fall (September, October, 
November). Silicate was relatively high toward the end 
of the spring season of upwelling, but decreased as the 
summer of 1963 approached and reached a minimum 
of about 3 microgram-atoms per liter in late July, 1963. 
Phosphate remained relatively stable with only minor 
increases in the spring. There was a consistent relation- 
ship between changes in the percentage saturation of 
oxygen, and pH. Nitrate showed an inverse relationship 
to both oxygen and pH, suggesting that changes in 
nitrate might be more closely related to biological 
processes than phosphate or silicate. In January, 1964, 


nitrate was minimal and was undetectable at the level 
of sensitivity of the analysis. 


In the waters along the shoreline of Southern Cali- 
fornia, light and carbon dioxide do not appear to limit 
the rate of reproduction of phytoplankton. Indeed, 
small intense local patches of phytoplankton have been 
observed during every month in which observations 
were made including the winter months of December 
and January when daylight is at a seasonal minimum. 
The pH range of 8.0-8.6, and the CO, range of 40-46 
ml/L recorded by Moberg, et. al. (81), indicate that 
carbon dioxide is always present in non-limiting quan- 
tities in the inshore waters. Temperature is not a 
critical factor in the growth of phytoplankton between 
5° and 25° C. as carbon fixation is principally a light- 
dependent reaction. Neither silicate-silicon or phos- 
phate-phosphorous have been determined to be present 
in any sample at less than 2 and 0.2 microgram-atoms 
per liter, respectively. Only nitrate has been found at 
essentially zero concentration. Data on absolute quan- 
tities of nutrients present in the water are perhaps 
suggestive, but the question of what constitutes a 
limiting quantity in the sea of any substance, either 
alone or in a synergistic system, is still an open question. 


The distribution of the surface sewage field is depicted 
in Figures 69 through 73 which show typical patterns 
of enrichment by phosphate and silicate of the area 
around the Orange County outfall. For this purpose 
phosphate and silicate were combined because their 
histories within the field during any one cruise were 
much the same. The effluent field generally flowed 
parallel to shore in an easterly direction but there was 
a seasonal variation. Flow to the northwest was ob- 
served more commonly in the winter, and flow to the 
east predominated in the summer. A similar directional 
pattern of inshore currents has been observed off 
Whites Point, Palos Verdes, some 30 miles northwest 
of the Orange County outfall (Malcolm Whitt, personal 
communication). 

The substances for which analyses were made and 
which were also contained within the sewage, may be 
divided into two general groups: 1) substances found 
in appreciable quantities in both the effluent wastes 
and in the receiving waters, e.g., silicate, phosphate 
and carbonate and 2) substances unique to the effluent 
waste discharge and found in negligible quantities in 
the receiving waters, e.g., ammonia and, experimentally, 
rhodamine-B. The difference between these two classes 
of substances is illustrated in Figure 74 where the 
relatively higher levels of SiO. and PO,-P are due to 
mixing with background waters which contain appreci- 
able quantities of those substances. However, Figure 
75 shows that the area receives a greater relative 
enrichment in ammonia than in either phosphorous or 
silicon. In this area, ammonia may offer a potentially 
useful enrichment to the receiving waters which are 
already sufficiently supplied with silicate and phosphate. 
As for the unknown and variable components of the 
effluent waste field, their possible role has yet to be 
investigated with respect to their biochemical effects. 

At this time, it appears that little nutritional benefit 
is gained from the addition of silicate and phosphate 
to waters where those nutrients may not limit the pro- 
ductivity of phytoplankton. Ammonia may provide 
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some useful nutriment as a supplement to the available 
fixed nitrogen normally found in the local receiving 
waters; however, the total contribution here is slight 
with respect to the total volume of inshore water along 
the coast as a whole. 

If effluent outfalls should contribute nutrients in 
sufficient quantity to influence the standing crops of 
phytoplankton, such an influence may be observed in 
several ways (Davis (82) ). The annual crop of phyto- 
plankton should steadily and significantly increase each 
year, or in the winter minima in phytoplankton should 
increase each year, or both may occur. Continuous 
records of the volumes of phytoplankton have been 
maintained for the last etght years by the Los Angeles 
City Bureau of Sanitation in connection with the 
routine monitoring of the Hyperion submarine outfall 
in Santa Monica Bay. These records were kindly made 
available for inspection by Mr. Norman Hume, Director 


of the Bureau of Sanitation. There has been no indica- 
tion of significant change; there has not been a signifi- 
cant increase in the standing crop of phytoplankton 
during the spring blooms, and the minima in succeeding 
winters fell to about the same values. This may be 
expected along an open coastline where a complete 
replacement of water may occur at least once each year 
at differing rates. Thus the oceanic situation may not be 
comparable to an essentially closed system such as 
Lake Erie (Davis (82) ) where the rates of eutrophica- 
tion exceeded those of turnover and exchange. 

This is but one of the many differences which exist 
between waste discharge into the marine environment 
and that found in a stream or lake. Furthermore, each 
marine environment will have its own unique features; 
the Pacific is generally rich in phosphorus and poor 
in nitrogen while the Mediterranean is poor in phosphor- 
us and rich in nitrogen. 


F. SUMMARY AND CONCLUSIONS 


1. DYE DIFFUSION STUDIES 


The principal objectives of the experiments were: 
(1) to verify some of the proposed multi-dimensional 
diffusion models as they pertain to the dispersion of 
dye from instantaneous and continuous point sources 
in the ocean, (2) to determine the effects of wind speed 
and water column stability on the parameters controll- 
ing the diffusion process, (3) to determine the applica- 
bility of the point source models in item (1) to the case 
of diffusion of a waste material from a large scale or 
volume source in the ocean, and (4) to propose, in the 
event the point source models did not apply, a solution 
for the problem of diffusion from a large scale instan- 
taneous or continuous source in the ocean. It was 
believed that the attainment of these objectives would 
aid considerably in an understanding of the basic 
parameters controlling the eddy diffusion process. In 
addition it was believed that information of this nature 
would aid considerably in the design of future waste 
outfall systems. 

A total of 32 diffusion experiments were conducted 
in near shore waters along the coast of Southern Cali- 
fornia. These experiments consisted of both instantane- 
ous and continuous release of tracer dye from a quasi- 
point source. In addition, one large scale experiment 
was accomplished in which an entire waste field was 
tagged with dye prior to discharge through a submarine 
outfall. The results of these experiments lead to the 
following conclusions: 


Relative Diffusion 


1. The three-dimensional point source models pro- 
posed by Gifford, which describe the space-time rela- 
tionship of mean and maximum concentration distribu- 
tion in a homogeneous and stationary turbulent field, 
were verified for the case of relative diffusion of dye 
patches at the sea surface. 

2. The rate of vertical diffusion, although much 
smaller than that in the corresponding horizontal 
directions, can contribute significantly to the overall 
diffusion process. The influence of vertical diffusion 
becomes increasingly important. on days of relatively 
high wind speeds (w + 8 knots) and/or low water 
column stability (S’ < 400). 

3. Horizontal diffusion models which neglect vertical 
diffusion can, for conditions of relatively low wind 
speed and/or high stability, be used to describe the 
decrease in maximum dye concentration with time. 
However, these models cannot be used to accurately 
predict the spatial distribution of dye concentration 
where the apparent rates of longitudinal and lateral 
diffusion differ considerably, as in the case of these 
experiments. 

4. The rate of vertical diffusion determined in these 
studies was demonstrated to be a function of average 
wind speed and average stability. At the same stability, 
the rate of vertical diffusion increased with increasing 
wind speeds. Conversely, at essentially the same wind 


speed, vertical diffusion proceeded at lower rates for 
high values of stability. 


5. The rate of longitudinal and lateral diffusion ap- 
peared to be influenced by wind speed but not by 
stability of the water column. The experiments con- 
ducted in the waste field indicated that the rate of 
longitudinal diffusion increased and lateral diffusion 
possibly decreased with increasing wind speeds. The 
latter trend was not well defined and additional data 
is required before a definite conclusion can be drawn 
in this respect. 


6. Although not specifically discussed in this report, 
it was not possible to establish a correlation between 
wind speed or individual diffusion parameters and mean 
current velocity. There was also no apparent correlation 
between the individual diffusion parameters, ¢,,y,2” 
themselves. 


Continuous Release Experiments 


1. The two-dimensional point source model proposed 
by Gifford which describes the maximum concentration- 
distance relationship along the center line of a steady 
plume in a homogeneous and stationary turbulent field, 
was found to be applicable to the results of the con- 
tinuous release experiments. 


2. Vertical diffusion, for conditions of relatively high 
wind speeds and/or low stability, contributed signifi- 
cantly to the diffusion and dilution of a steady dye 
plume. 


3. The one-dimensional (horizontal) diffusion models 
proposed by Batchelor and by Schonfeld were found 
to be applicable to the experiments where the rate of 
vertical diffusion was suppressed by high water column 
stability and/or low wind speed. 


4. The results of the single full scale experiment in 
which the entire waste field was tagged, showed that 
diffusion models of a point source type, based on an 
assumed Gaussian concentration distribution, cannot 
accurately describe the observed diffusion phenomenon. 
The principal reasons for the discrepancy are: (a) the 
lateral dye concentration distribution within the surface 
field at Orange County is non-Gaussian, and (b) the 
initial size of the waste field in no way resembles a true 
point source. 

5. Equations for one- and two-dimensional diffusion 
from instantaneous and continuous volume sources are 
presented as a realistic model of the observed large 
scale diffusion phenomenon. These equations are based 
on the assumptions of an initial Gaussian distribution 
of material in a homogeneous and stationary turbulent 
field, as originally proposed by Gifford. The proposed 
equation for diffusion from a continuous source was 
derived from Gifford’s instantaneous model by analogy 
with similar diffusion models. The assumption of a 
Gaussian material distribution appears to be reasonable 
for outfalls incorporating a more efficient diffuser sys- 
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tem than that presently found at the Orange County 
site. 


General 


1. The data from both instantaneous and continuous 
release experiments conducted within a waste field 
indicate that the spatial distribution of dye concen- 
tration at any time (or distance) is approximately 
Gaussian. This is especially true at small values of 
dispersion time (< 30 minutes). In some experiments 
the observed dye distributions were skewed to the right 
or left of the maximum concentration; however, it 
appears the Gaussian distribution can be used in most 
cases as a good working rule. 

2. The data from these experiments indicates that 
Batchelor’s dimensional predictions regarding the rela- 
tionship between lateral variance and diffusion time 
(or distance) could not be verified. This implies that, 
providing the rate of energy dissipation was approxi- 
mately constant, the “4/3 law” relating the lateral 
coefficient of eddy diffusion as a function of average 
eddy scale, does not hold in the particular oceanic areas 
studied in these experiments. 


2. PRIMARY PHYTOPLANKTON PRODUCTION 
AND PRODUCTIVITY 


During this project the production and productivity 
of the phytoplankton has been used to measure the 
biological effect of sewage disposal in the open coastal 
marine environment. Surveys of an area affected by 
sewage showed a pattern of low surface productivity 
near the point of discharge. An increase in productivity 
to values well above background was found downstream 
after a time of about 6 hours. Maximum values occurred 
at about 10 hours. This was followed by a decrease 
towards background levels. The spatial location of these 
features depended on the direction and magnitude of 
the current and the rates of dilution. 

A seasonal pattern of high productivity in spring 
and summer (9.9 metric tons of C fixed per day in an 
affected area of 60.2 Km?) and low productivity in fall 
and winter (4.6 metric tons) was found. Red tide con- 
ditions resulted in increased productivity to as much 
as five times the normal seasonal values. 

Surveys along a roughly similar coastal area that was 
presumed to be free of sewage showed only the expected 
pattern of higher production inshore. There was evi- 
dence of a seasonal pattern similar to that in the area 
affected by waste. Surface productivity in the unaffected 
area was roughly half of that in the area containing 
waste. 


The same sequence in the pattern of productivity 
found in the affected area was produced by experimental 
exposure of subsamples of a uniform phytoplankton 
population to serially diluted effluent. After fourteen 
hours of storage no material changes were observed in 
the rates of photosynthetic uptake in individual sam- 
ples. The effect of the sewage therefore appears to be 
closely related to dilution. 


The relationship of productivity to specific features 
of the environment in the sewage field was explored. 
Although the sewage itself varies in content and quan- 
tity, the general pattern of productivity in relationship 
to the outfall has been consistent. Therefore, the 
attempt to relate the pattern to specific limiting or 
inhibiting factors was directed towards those compon- 
ents of the sewage that appeared to be always present 
and for which convenient analytical techniques were 
available. Among the conditions imposed on the re- 
ceiving waters are the reduction of water transparency 
which limited subsurface production near the boil, and 
the added load of nutrients which apparently are 
neither limiting nor stimulating since they showed no 
consistent statistical relationship to productivity. There 
was no evidence that the nutrients were limiting in the 
waters free of sewage except in the areas greater than 
3 miles from shore. 


3. ULTRAVIOLET LIGHT AS A MEANS OF 
TRACKING A SEWAGE FIELD 


It has been shown that the ultraviolet absorption of 
effluent at 230 millimicrons from the Orange County 
outfall was considerably higher than that of background 
waters. The distribution over the waste field was similar 
to that of ammonia, although this does not imply any 
direct relationship between the two except that they 
are both measures of physical dilution. 

Measurement of UV absorption is simple and rapid 
with the Beckman Model DB spectrophotometer. 


4. CHEMICAL PROCEDURES 


For use in conjunction with both dilution and pro- 
ductivity studies, analyses were made for salinity, 
dissolved oxygen, pH and alkalinity, phosphate, silicate, 
nitrate, nitrite and ammonia. Of these, ammonia proved 
to be the most useful of the chemical tags with which to 
describe the extent and dilution of the effluent field 
although the analytical procedure leaves much to be 
desired for routine or field use. Much remains to be 
done in the development of suitable techniques for 
most of the substances listed. 
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LIST OF SYMBOLS 


Definition 


Major and minor semi-axes of a semi-ellipsoid. 

Variance of tracer concentration at the point 
of origin. | 

Mean tracer concentration as a function of time 
(or distance) and location within a waste field. 

Initial maximum tracer concentration at the 
point of origin. 

Maximum tracer concentration as a function of 
time or distance. 

Average depth of tracer or waste field. 

Minimum dilution of tracer as a function of time 
or distance. 

Variance of the center of mass of tracer sub- 
stance. 

Water column stability. 

Approximate water column stability. 

Rate of energy dissipation per unit mass. 

Variance of diffusing particles in given co- 
ordinate direction. 

Gravitational acceleration. 

Constant. 


/ 
xy Apparent coefficient of eddy diffusion in co- 


ordinate directions. 
Rate constant. 


+*’Y Lagrangian time scale of turbulence in coor- 


dinate directions. 

Average eddy scale. 

Total amount of tracer substance released from 
an instantaneous source. 

Total amount of tracer as computed from mass 
balance analysis. 


Symbol 
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X,Y,Z 


(115 ) 


Definition 


Total amount of tracer per foot of depth dis- 
charged from an instantaneous source. 

An exponent. 

In situ water temperature. 

Most probable diffusion velocity. 

Rate of tracer discharged from a continuous 
source. 

Average rate of tracer discharged from a con- 
tinuous source. 

Average displacement of tracer particles with 
respect to center of tracer mass. 

Initial particle separation. 

Density. 

Bex pl 0%, 

Standard deviation of particle or tracer concen- 
tration distribution as a function of time or 
distance. 

Average variance of particle or tracer concen- 
tration distribution as a function of time or 
distance. 

Density factor. 

Diffusion time. 

Observed diffusion time. 

Average speed of surface current. 

Average wind speed. 

Weighted mean diffusion velocity. 

Coordinate axes with x in the direction of the 
mean current, y perpendicular to x in the 
lateral direction, and z in the vertical direc- 
tion. 
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